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Abstract. The covariant density functional theory with a few number of parameters 
has been widely used to describe the ground-state and excited-state properties for 
the nuclei all over the nuclear chart. In order to describe exotic properties of 
unstable nuclei, the contribution of the continnum and its coupling with bound states 
should be treated properly. In this Topical Review, the development of the covariant 
density functional theory in continuum will be introduced, including the relativistic 
continuum Hartree-Bogoliubov theory, the relativistic Hartree-Fock-Bogoliubov theory 
in continnum, and the deformed relativistic Hartree-Bogoliubov theory in continuum. 
Then the descriptions and predictions of the neutron halo phenomena in both spherical 
and deformed nnclei will be reviewed. The diffuseness of the nuclear potentials, nnclear 
shapes and density distributions, and the impact of the pairing correlations on nuclear 
size will be discussed. 
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1. Introduction 

The building blocks of nucleus are protons and neutrons. An arbitrary combination of 
protons and neutrons doesn’t necessarily make a nucleus. Among thousands of nuclides 
expected to exist, only 288 of them are stable or at least with half-lives longer than the 
expected life of the Solar System, which form the valley of stability in the nuclear chart. 
These nuclides are sandwiched by about 3,000 nuclides produced in laboratories [1]. 
The nuclear density functional theory (DFT) has predicted that about 7,000 nuclides are 
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bound with respect to neutron or proton emission [2, 3, 4], If the continuum contribution 
is included, as demonstrated in Ref. [5] from O to Ti isotopes, this number may exceed 
10,000. In addition, there are also nuclides (resonances) which, though beyond the drip 
lines, can have lifetimes long enough to be measured and investigated. 

The nuclides lying above the valley of stability in the nuclear chart are proton-rich 
and the proton drip line has been identified experimentally up to protactinium [Z = 91). 
However, for those lying below the valley of stability, which are neutron-rich, the border 
of the nuclear territory is known only up to oxygen [Z = 8). Experimental exploration of 
very neutron-rich nuclei is extremely challenging because of very low production rates 
in studies involving the fragmentation of stable nuclei, as well as the separation and 
identification of the products. 

In nuclear physics, the study of the properties of exotic nuclei — nuclei with extreme 
numbers of proton or neutron — is one of the top priorities as it can lead to new insights 
on the origin of chemical elements in stars and star explosions. Although the radioactive 
ion beams (RIB) have extended our knowledge of nuclear physics from stable nuclei to 
exotic ones far away from the valley of stability, it is still a dream to reach the neutron 
drip line up to mass number A ~ 100 with the new generation RIB facilities developed 
around the world, including the RIKEN Radioactive Ion Beam Factory (RIBF) in 
Japan [6], Facility for Antiproton and Ion Research (FAIR) in Germany [7], Second 
Generation System On-Line Production of Radioactive Ions (SPIRAL2) at GANIL in 
France [8], Facility for Rare Isotope Beams (FRIB) in USA [9], Rare Isotope Science 
Project (RISP) in Korea [10], Gooler Storage Ring (GSR) at Heavy Ion Research Facility 
in Lanzhou (HIRFL) in Ghina [11, 12], etc. 

The development of RIB facilities around the world significantly stimulates the 
study of exotic nuclei far from the valley of stability [13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23]. New and exotic phenomena were observed in nuclei close to drip lines such 
as neutron or proton halos [24], changes of nuclear magic numbers [25], the island of 
inversion [26], pygmy resonances [27], etc. More exotic nuclear phenomena have also 
been predicted, e.g., giant halos [28, 29], shape decoupling between core and halo [30, 31], 
etc. 

In halo nuclei, the extremely weak binding of valence nucleons leads to many new 
features, such as the coupling between bound states and the continuum due to pairing 
correlations and the very extended spatial density distributions. Therefore, one must 
consider properly the asymptotic behavior of nuclear densities at a large distance from 
the center and treat in a self-consistent way the discrete bound states, the continuum, 
and the coupling between them in order to give a proper theoretical description of the 
halo phenomenon [20, 32]. There are mainly two key points to achieve a self-consistent 
description of halo nuclei with density functional theories: the pairing correlations and 
the suitable basis. 

In the mean field descriptions of nuclear halos, nuclear pairing plays an important 
role [33, 34, 35, 28, 29]. The conventional BGS method for the pairing is not appropriate 
for exotic nuclei. To include properly the contribution of continuum states, the 
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Bogoliubov transformation has been justified to be very useful [36, 37]. By solving 
the nonrelativistic Hartree-Fock-Bogoliubov (HFB) [36, 37] or the relativistic Hartree- 
Bogoliubov (RHB) [34, 38, 39] equations in coordinate (r) space, one can fully take into 
account the coupling to the continuum. As an alternative, one may hrst locate resonant 
states in the continuum and compute the resonance parameters with various methods 
[40, 41, 42, 43, 44, 45, 46], then by using the resonant-BCS (rBCS) approach [47], include 
the contribution of these resonant states and study halo phenomena [48, 49, 50]. The 
Green’s function method [51], which can take into account properly the asymptotic 
behavior of continuum states, has also been adopted for solving the HFB [52, 53, 54], 
the self-consistent Skyrme-HFB [55, 56, 57] or RMF [58] equations in the coordinate 
space. 

The solution of the coupled differential equations of HFB and RHB theories is 
relatively easier for spherical systems with local potentials, where one-dimensional 
Numerov or Runge-Kutta methods can be applied. This is also true even for nonlocal 
problems where the hnite element method (FEM) can be used [39] . A more widely used 
method for solving such equations is to expand the single particle wave functions in a 
basis and this is more convenient particularly for deformed systems. In the past decades, 
the harmonic oscillator (HO) basis has been used with a great success for stable nuclei 
in both nonrelativistic and relativistic mean held approximations such as the Skyrme 
Hartree-Fock, Hartree-Fock-Bogoliubov, relativistic Hartree, and relativistic Hartree- 
Bogoliubov theories for deformed or nonlocal systems [59, 60, 61]. However, due to the 
incorrect asymptotic behavior of the HO wave functions, the expansion in a conventional 
HO basis is not proper for the description of exotic nuclei [62]. Considerable efforts 
have been made to develop mean held models in a basis with an improved asymptotic 
behavior at large distances, e.g., a transformed HO basis [63, 64] or a Woods-Saxon 
(WS) basis [65]. 

The DFT with a small number of parameters is widely used to study the ground- 
state and excited-state properties of the nuclei all over the nuclear chart. In particular, 
the covariant density functional theory (CDFT), incorporating the Lorentz symmetry 
in a self-consistent way, has received wide attention [66, 67, 68, 20, 69, 70]. 

There exist a number of attractive features to describe atomic nuclei in the 
relativistic framework. The most obvious one is the natural inclusion of the nucleon 
spin degree of freedom, resulting in the nuclear spin-orbit potential that emerges 
automatically with the empirical strength in a covariant way. Moreover, it provides 
a new saturation mechanism for nuclear matter [71], reproduces well the measured 
isotopic shifts in the Pb region [72] , reveals more naturally the origin of the pseudospin 
symmetry [73, 74] as a relativistic symmetry [75, 76, 77, 78, 79, 80, 81, 82], and predicts 
the spin symmetry in the anti-nucleon spectrum [83, 84]. It can also include the nuclear 
magnetism [85], i.e., a consistent description of currents and time-odd helds, which plays 
an important role in the nuclear magnetic moments [86, 87, 88, 89] and nuclear rotations 
[90, 91, 92, 93, 70], etc. 

The CDFT is a reliable and useful model for nuclear structure study in the whole 
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nuclear chart, and therefore, it is natural and necessary to describe the halo nuclei based 
on the CDFT. 

Since the discovery of the neutron halo, i.e., the anomalously large nuclear radius in 
^^Li [24] which deviates from the conventional law on the mass dependence of the 
nuclear radius, considerable efforts have been undertaken to understand this interesting 
phenomenon in exotic nuclei close to the neutron drip line. 

The relativistic mean field models have been used to study nuclear halo phenomema 
since early 1990s [94, 95, 96, 97, 98, 99, 100]. By incorporating the Bogoliubov 
transformation in the relativistic Hartree theory, the relativistic continuum Hatree- 
Bogoliubov (RCHB) theory has been developed and extensively used to study the 
halo phenomeon in spherical nuclei [34, 38, 28, 101, 39]. The RCHB theory was also 
extended to including the hyperon and used to study neutron halo in A-hypernuclei 
[102, 103, 104, 105, 106]. The relativistic Hartree-Fock-Bogoliubov theory in continuum 
(RHFBc) theory has been developed and used to study the influence of exchange 
terms on the formation of neutron halos [107]. Over the past years, lots of efforts 
have been made to develop a deformed relativistic Hartree-Bogoliubov theory in 
continuum (DRHBc) [108]. Halo phenomena in deformed Ne and Mg isotopes have 
been investigated with the DRHBc theory [30, 31, 109]. 

In this Topical Review, recent progresses on the application of the CDFT for the 
neutron halo phenomena in spherical and deformed nuclei will be presented. As an 
extensive review is available [20], the RCHB theory and its application will be briefly 
mentioned here for completeness. In Section 2, we give the formalism for CDFT in 
continuum, including the general framework of CDFT, the formalism for RHB theory, 
and the RCHB, RHFBc and DRHBc theories. In Section 3, the applications of the 
CDFT for halos and giant halos in spherical nuclei will be reviewd. We will discuss the 
RCHB description of the neutron halo in ^^Li, the prediction of giant halos in Zr and Ca 
isotopes, the continuum contribution and the extension of nuclear landscape, halos from 
relativistic and non-relativistic approaches, and halos with RHFBc theory. The densities 
and potentials in exotic nuclei become more diffuse and the surface diffuseness has 
important influences on spin-orbit splittings, single particle level structure, and nuclear 
masses. Recent progresses on these topics will be presented in Section 4. In Section 5, 
we will discuss the influence of pairing correlations on the nuclear size. The necessity to 
carry out a self-consistent calculation will be shown in order to understand the impact 
of the pairing correlations on the nuclear size. In Section 6, the study of neutron halo in 
deformed nuclei will be reviewed, including the neutron separation energies, quadrupole 
shapes, density distributions, and radii. The generic conditions for the occurrence of 
halos in deformed nuclei and the predicted shape decoupling between the core and the 
halo will be discussed. Finally we give a summary and discuss perspectives in Section 7. 
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2. Covariant density functional theory in continuum 


2.1. General framework of CDFT 


The covariant density functionals represent a number of relativistic mean field (RMF) 
and relativistic Hartree-Fock (RHF) models based on the framework of quantum 
hadrodynamics (QHD) [66]. The advantages in using covariant functionals include the 
natural inclusion of the nucleon spin degree of freedom, the consistent treatment of 
isoscalar Lorentz scalar and vector self-energies which provides a unique parametrization 
of time-odd components of the nuclear mean held, a natural explanation of the empirical 
pseudospin symmetry, and a new saturation mechanism of nuclear matter with a 
distinction between scalar and four-vector nucleon self-energies [67, 110, 111, 68, 20, 
112, 82], 

The conventional CDFTs are based on the hnite-range meson-exchange 
representation, in which a nucleus is described as a self-bound system of Dirac nucleons 
interacting with each other via the exchange of mesons through an effective Lagrangian 
density, which is associated with the nucleon (-0), the a-, u-, p-, and vr-meson helds, and 
the photon held (A), 


C — C]\f -|- Cfj + + Cp + + Ca + -^ 7 , ( 1 ) 

containing the Lagrangian densities for the nucleon 

C.N = 'ip - M) fj, ( 2 ) 


for the cr-, cn-, p-, and 7r-mesons, and the photon A, 
jC.a= + ]^d^adpa - 

^P= - + ^rnlp^ ■ Pp, 

R-k= + ^dpTl ■ <9^7? - ^mlTi ■ 7f, 

Ra= - 

and for the interactions between nucleon and meson helds. 


Cl 




f 

-9a(^ - gcoYujp - gpYr ■ pp - ^ 


-—Ibl^dpTi -r-ey^ ^ Ci.p 
m.,r 2 


fj. 


( 3 ) 


In above Lagrangian densities, M denotes the mass of nucleon, and m^j {go), 
{gui), (5'p) fp)i and (/^r) are respectively the masses (coupling constants) for a-, ca-, 
p-, and TT-mesons. The p-tensor coupling term, which is practically negligible at Hartree 
approximation, can signihcantly improves the descriptions of nuclear shell structures at 
Hartree-Fock approximation [113]. 
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The field tensors of the vector mesons and R and the electromagnetic field F are 
defined as, 

(4) 


The isoscalar-scalar cr-meson provides the long-range attractive part of the nuclear 
interaction whereas the short-range repulsion is governed by the isoscalar-vector u- 
meson. The photon field A^^{x) accounts for the Coulomb interaction and the isospin 
dependence of the nuclear force is described by the isovector-vector p-meson. The vr- 
meson field and the rho-tensor coupling contribute only if the Fock terms are included. 

In the following, the vectors in the isospin space are denoted by arrows and the 
space vectors by bold type. Greek indices p and i' run over the Minkowski indices 0, 1, 
2, 3 or f, X, y, z, while Roman indices i, j, etc. denote the spatial components. 

The variation of the Lagrangian L = J (RxC{x) with respect to nucleon (-^) and 
meson fields (0j) leads to the Dirac equation of nucleon, 

+ M + S) ij{x) = 0, (5) 


with the self-energy S to be self-consistently determined, and the Klein-Gordon 
equations for the meson fields (0 = a,and if), 

+ ml)(t) = ( 6 ) 


with the source terms 
' 

S^= < 


0 = a, 

Qpi’lt.T'ilJ + (!> = Pp, 

0 = TT, 


( 7 ) 


and the Proca equation for the photon field (R), 

d''F^^ = ( 8 ) 

Under the no-sea and mean field approximations, the Dirac Hartree-Fock equation 
can be derived as 


J dr'h{r,r')'ijj{r') = 


( 9 ) 


where e is the single particle energy (including the rest mass) and the single particle 
Dirac Hamiltonian h{r, r') contains the kinetic energy the direct local potential 
hP, and the exchange non-local potential hP, 


pkin(^, r') = [a-p + /3M] (5(r - r'), 
h°(r, r') = [ST(r )75 + So(r) + /3S5(r)] 5{r - r'), 
YG{r,r') Yp^r.r') 

Xcir^r’) XF{r,r') 




( 10 ) 

( 11 ) 

( 12 ) 
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In the above expressions, the local scalar S 5 , vector Eg, and tensor self-energies 
contain the contribntions from the direct (Hartree) terms and the rearrangement terms, 
while the nonlocal self-energies Xg-, Xp, Yq, and Yp come from the exchange (Fock) 
terms. More details can be fonnd in Ref. [114]. 

If one sticks to the Hartree level (RMF), for a nucleus with time reversal symmetry, 
the Dirac equation (5) is reduced as 

h'ip{r) = e'lpkir), (13) 

where 

h = cx-p + l3[M + S 5 (r)] + So(r), (14) 

is the single particle Hamiltonian. 

More recently, the RMF framework has been reinterpreted by the relativistic Kohn- 
Sham density functional theory, and the functionals have been developed based on the 
zero-range point-coupling interaction [115, 116, 117, 118]. In fact, by solving formally 
the Klein-Gordon equations for mesons and taking the large meson mass limit, one can 
derive point-coupling functionals in which the meson in each channel (scalar-isoscalar, 
vector-isoscalar, scalar-isovector, and vector-isovector) is replaced by the corresponding 
local four-point (contact) interaction between nucleons. 

The basic building blocks of CDFT with point-coupling vertices are 

(i/^OTV'), C>G{l,r}, F e { 1 , 7 ^, 75 , 757 ^, CT^^}, (15) 

where ijj is the Dirac spinor held of nucleon, r is the isospin Pauli matrix, and F 
generally denotes the 4x4 Dirac matrices. There are ten of such building blocks 
characterized by their transformation characteristics in isospin and Minkowski space. 
A general effective Lagrangian can be written as a power series in and their 

derivatives, with higher-order terms representing in-medium many-body correlations. 
For details, see Refs. [112, 115, 116, 117, 118]. 

The point-coupling model has attracted more and more attention owing to the 
following advantages. First, it avoids the possible physical constrains introduced by the 
explicit usage of the Klein-Gordon equation to describe the mean meson helds, especially 
the hctitious a meson. Second, it is possible to study the role of naturalness [119, 120] 
in effective theories for nuclear structure related problems. Third, it is relatively easy 
to include the Fock terms [121], and provides more opportunities to investigate its 
relationship to the nonrelativistic approaches [ 122 ]. 

2.2. Relativistic Hartree-Fock-Bogoliubov theory 

Pairing correlations are crucial in the description of open-shell nuclei. For exotic nuclei, 
the conventional BGS approach turns out to be only a poor approximation [37, 20]. 
Starting from the RMF Lagrangian density, a relativistic theory of pairing correlations 
in nuclei has been developed by Kucharek and Ring [123, 67], 



( 16 ) 
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with 

( h{r,r') - X5{r-r’) A(r,r') \ 

H{r,r) = , (17) 

y A(r,r') —h{r,r') + X6{r — r') J 

where the chemical potential A is introduced to preserve the particle number on the 
average. 

The single particle Hamiltonian h{r, r') is derived from Eq. (5) with the 
retardation effects neglected, e.g., the Hamiltonian given in Eq. (13) under the Hartree 
approximation. The pairing potential can be written as, 

Aa(r, r') = -\Y^ (r, r') fi:^(r, r'), (18) 

^ /3 

where the pairing tensor k is 

^^a{r,r') ='4)vSr)*'4)uSr')- (19) 

The pairing force is either taken as a density-dependent two-body force in a zero range 
limit, 

V(r, r') = VoS{r - r')^ (1 - cr ■ cr') ^1 - , (20) 

with an adjusted strength Vq, or as the pairing part of the Gogny force [124], 

V{r, r) = (H7 + BiP^^ - HiP^ - M^P^P^) , (21) 

1 = 1,2 

with the parameters /ij, H7, Hj, Hi and Mj. It has been shown that after proper 
renormalization, the above density-dependent force of zero range and the finite range 
Gogny force produce essentially the same results [35]. 


2.3. Relativistic continuum Hartree-Bogoliubov theory 

If the Fock terms are neglected, as it is usually done in the covariant density functional 
theory, the Dirac Hartree-Bogoliubov equation for the nucleons still takes the same form 
as that given in Eq. (16). For spherical nuclei, the quasiparticle wave functions read, 

, ^ 1 / iGUr)yU0.ct>)\ 

, ^ 1 / \G\,{r)y]^{e,ct>)\ 

The angular part of the RHB equation can be integrated out and the radial RHB 
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equations are derived as 

+ -Gu{r) -{E + \- So(r) + Ss(r))Ff/(r) +r [ r'dr'Apir, r')Fv{r') = 0, 

dr r J 

- -Fu{r) + {E + \- So(r) - S 5 (r))Gc/(r) + r [ Edr'Acir, r')Gy (r') = 0, 
< ^r( \ ^ r (2^ 

—5^ + -Gv{r) + {E-X + So(r) - Es{r))Fv{r) + r / r'dr'Apir, r')Fu{r') = 0, 

dr r J 

dFv{r) _ So(r) + Sg(r))Gv(r) + r f r'dr'Acir, r')Gu{r') = 0. 

^ ctr r J 

The pairing potentials are calculated as 


AG(r, /) = 15^ ri7"(r,r')9,. [GprjGpr') + Gpr')G‘(r)], (25) 

AF(r,rO = i^nt’“(nOfe ^[f5(r)fv(’'') + fi(Of^(r)|. (26) 

k' £fe >0 

If the zero-range pairing force is used, the above coupled integro-differential equations 
are reduced to differential ones, which can be directly solved in coordinate space [125]. 


2.4- Relativistic Hartree-Fock-Bogoliubov theory in continuum 


For spherical systems, the solution of the RHFB equations, i.e., the Dirac spinor 'ipUa 
and 'ipVa) can be written similarly to Eqs. (22) and (23), 


^ ^ ^ I '^GuSr)yf^^^{r) 

" V -FuSr)y'f^mS^) 


Mr) = 




Jama ' 


(27) 


(28) 


The RHFB equations (16) are then reduced to the coupled integro-differential equations, 

Gi,.(r)-(B„+A-E_)F(,.(r) 


|- + ^ + E^ 

dr r 


+ XuA'r) + r J r'<ir'A„(r,r')JV(r') = 0, 
FG.(r) + (B„ + A-E+)GG.(r) 

Jrdr'Aair,r')Gv^{r’) ^ 0, 

+ Et GK,(r) + (E„-A + E_)F,nr) 


(29) 


d „ 

-- 

dr r 


d K, 


dr r 

+ XvSr) + r j r'dr'Aair,r')Fu,{r') = 0, 
d Ka 


dr r 


n.(r)-(F„-A + E+)GK,(r) 


n,(r) + r / r'dr'A^(r,r')G,j,(r') = 0, 
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where Ea are the quasi-particle energies (without the rest mass), and the local self¬ 
energies S_|_ and S_ are 


S_|_ = So -|- S 5 , (30) 

S_ = So - S 5 - 2M. (31) 


More details can be found in Ref. [114], 

The pairing potentials A(j(r, r') in Eq. (29) can be expressed as 

Aa(r,r') = -^1/J(r,r')fi:fe(r,r'), (32) 

b 

where the pairing tensor K{r, r') reads 


^^b{r,r') = [Gu,{r)Gv,{r') + Fu,{r)Fv,{r')\ + ]^jl [Gv,{r)Gu,{r') + Ei4(r)Ff/^(r')], (33) 


where = 2j -|- 1. Details of the pairing interaction matrix element can be found 
in Ref. [39]. 

In contrast to the RHB approach with (5-force in the pairing channel where the 
radial equations (29) become differential equations, in RHFB theory the radial equations 
are fully integro-differential. The integral terms may arise from the Fock terms or the 
pairing channel if the finite-range pairing force is used. In coordinate space, it is difficult 
to solve such equations by the localization procedure adopted in Refs. [126, 127] when 
solving the relativistic Hartree-Fock equations. It is more convenient to solve them by an 
expansion of the Dirac-Bogoliubov spinors in an appropriate basis. This has been done 
by using the Dirac Woods-Saxon (DWS) basis [65]. This basis has been constructed for 
the investigation of weakly-bound nuclei. The set of DWS basis functions 


^ 0}, (34) 

are eigenfunctions (with eigenvalues Sh) of a Dirac equation with Woods-Saxon-like 
potentials for So(r) ±S 5 (r). They are determined by the shooting method in coordinate 
space within a spherical box of size i?max [128]. 

The U and V components of the Dirac-Bogoliubov spinors (27) can be expanded 
as, 

Np No 

^ UpQp + ^ Udgd, (35) 

p=l d=l 

Np Nd 

p=l d=l 

where Np and Np respectively correspond to the numbers of positive [Sp > 0) and 
negative [sd < 0) energy states in the DWS basis. Obviously, because of spherical 
symmetry the quantum number k is preserved, i.e., the RHFB equations can be solved 
for each value of k and the sums in the expansion (35) run only over states with the 
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same k. For a fixed value of k we have the radial basis spinors 


9p(r) 

9d{r) 


Gp{r) \ 

Fr(r) J 

Gd(r) A 

FdM J 


(36) 

(37) 


where the sub-indices p and d correspond to the number of nodes of the basis functions 
Gp for positive energy and for negative energy. 

In the DWS basis (35) the radial RHFB equations (29) are transformed to a matrix 
eigenvalue problem, 


^,-A A 

A -H + \ / \ V / \V J 


(38) 


where H and A are {Np + No) x {Np + iV/))-dimensional matrices, U and V are the 
column vectors with Np -h No elements. From the expressions of the single particle 
Hamiltonian h and pairing potential A given in the previous part we obtain the matrix 
elements of H and A as 


Hnn' = j drGnir) F^,{r) + j drFn{r) G'„K’^X39) 

, = J dr [Gnir)GAr)J:pir) + F4r)F„,(r)S_(r)] + 



(40) 

(41) 


Ann' = J dr j dr'A,,{r,r')[Gnir)Gn'ir') + Fnir)Fn>ir')], (42) 

where n and n' run over the radial quantum numbers of the DWS basis states in Eq. (37) 
with both positive energies {n,n' = p) and negative energies {n,N = d). 


2.5. Deformed relativistic Hartree-Bogoliubov theory in continuum 

For deformed nuclei the solution of HFB or RHB equations in r space is numerically 
a very demanding task. Considerable efforts have been made to develop mean held 
models either in r space or in a basis with an improved asymptotic behavior at 
large distances [129, 63, 130, 65, 131, 132, 133]. In particular, an expansion in a 
Woods-Saxon (WS) basis was proposed as a reconciler between the HO basis and the 
integration in coordinate space [65]. Woods-Saxon wave functions have a more realistic 
asymptotic behavior at large r than the harmonic oscillator wave functions. A discrete 
set of Woods-Saxon wave functions is obtained by using box boundary conditions to 
discretize the continuum. It has been shown in Ref. [65] for spherical systems that the 
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solution of the relativistic Hartree equations in a Woods-Saxon basis is equivalent to 
the solution in coordinate space. The Woods-Saxon basis has been used recently for 
both nonrelativistic [134] and relativistic Hartree-Fock-Bogoliubov theories [114] with a 
hnite range pairing force. 

Over the past years, lots of efforts have been made to develop a deformed relativistic 
Hartree (RH) theory [135] and a deformed relativistic Hartree-Bogoliubov theory in 
continuum (DRHBc theory) [108]. As a hrst application, halo phenomena in deformed 
nuclei have been investigated within the DRHBc theory [30]. The detailed theoretical 
framework is available in Ref. [31]. 

For axially deformed nuclei with spatial reflection symmetry, the potentials T,s{r) 
and Flo(r) in Eq. (14) and various densities can be expanded in terms of the Legendre 
polynomials [136], 

fir) = '^fxir)Px{cos9), A = 0,2,4, •••, (43) 

A 

with 

2A -1-1 

fxir) = —^— J d{cose)f{r)Px{cose). (44) 

The quasiparticle wave functions 'ipUk 'ipVk are expanded in terms of spherical 
Dirac spinors (pn.Kmi'f'sp) with the eigenvalues obtained from the solution of a 
Dirac equation containing spherical potentials and l/^°^(r) of Woods-Saxon 

shape [65, 128]: 


f^Uuirsp)= ^u^^l^^TnKmirsp), 

nn 

(45) 

i^v^rsp) = 

(46) 

The basis wave function reads 


. . 1 ( *G„«(r)3^j^(Ds) \ 

TnumyS) — 1 77 / /'O ^ ) ’ 

r \ -Fn^ir)y)^{ns) J 

(47) 


where GnKif')/r and FnKi'r)lr respectively the radial wave functions for the upper 
and lower components. The state ^nK.mi'f'sp) is the time reversal state of (pnumi'f'sp). 
The spherical spinor (pnKm is characterized by the radial quantum number n, angular 
momentum j, and the parity tt; the latter two are combined to the relativistic quantum 
number k = 7 r(—l)N^/^(^j -|- 1/2) which runs over positive and negative integers 
K = ±1,±2,---. For the spinor spherical harmonics and ^ = J + |sign(fi;) 
and I = j — |sign(fi;). These states form a complete spherical and discrete basis in Dirac 
space (see Ref. [65] for details). Because of the axial symmetry the 2 ;-component m of 
the angular momentum j is a conserved quantum number and the RHB Hamiltonian 
can be decomposed into blocks characterized by m and parity tt. For each (m7r)-block, 
solving the RHB equation is equivalent to the diagonalization of the matrix 


>f-A B \(ldk\ 
e* -A’ + \)\Vk) 



( 48 ) 
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(mir) 

^k,{nK) 


*K)(«'»')) = ((uftmirl/iDln'ft'mTT)), 
= ((uftmir|A|n'K'm7r)). 
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(49) 

(50) 

(51) 


Further details are given in Ref. [31]. 

After solving the RHB or RHFB equations self-consistently, one gets the 
qnasiparticle wave functions and meson helds of a nucleus. Then various physical 
quantities can be calculated. Here we take the deformed RHBc as an example and show 
how the total energy, the radius, and the deformation parameters are calculated [31]. 

The total energy of a nucleus is 


where 


E — F'nucleon + A’fr + A’u + F'p + F'c + A'c.m. 

~ 2(A — F'fc)'yfc^(m7r) ~ -^pair 

d^^[9a(^ir)ps{r) + U{a)] 

- \ J d^rgpp{r) [p^{r) - p^ (r)] 

- ^ / d^rAop^{r) + 

= J d^rV^{r)V,,{r) = ’ 


(52) 

(53) 


and U{a) is the nonlinear self-coupling term of a meson. 

For a zero-range force the pairing held A(r) is local, and the pairing energy is 
calculated as 


h/pair 


1 

2 


/ 


d^rK,{r)A{r). 


(54) 


Note that in the DRHBc theory, a smooth cutoff with two parameters, similar to the soft 
cutoff proposed in Ref. [137], has been introduced to regularize the zero-range pairing 
force; for the details, the readers are referred to Refs. [30, 31]. The center of mass 
correction energy 


Er.m. 


2Am 


A-), 


(55) 


is calculated after variation with the wave functions of the self-consistent solution [138, 
139] or in the oscillator approximation 


Bc.m, = X 41 X MeV, 


(56) 
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The root mean square (rms) radius is calculated as 


Rt 


= (r2)V2 = 


d?r 


1/2 


dr [r^pl,x=o{r)\ 


1/2 


(57) 


where r represents the proton, the neutron, or the nucleon. The rms charge radius is 
calculated simply as = r‘^ + 0.64 fm^. The intrinsic multipole moment is calculated 
by 


Qt,X 


IQtt 
2A + 1 
Stt / 


{r^Yxo{9,4>)) = 2{r^Pxm 
dr [r^pl,xi.r)] ■ 


2X + 1 

The quadrupole deformation parameter is obtained from the quadrupole moment by 

V^Qt,2 


(58) 


3Af,{r2) ’ 

where Nr refers to the number of neutron, proton, or nucleon. 


(59) 


3. Neutron halos and giant halos 

Halo is one of the most interesting exotic nuclear phenomena. There are many new 
features in the halo nuclei because of their large spatial distribution. To give a proper 
theoretical description of the halo phenomena, it is crucial to treat properly the coupling 
between bound states and the continuum due to pairing correlations by either the hnite 
range or the zero range pairing forces and the very extended spatial density distributions. 
By assuming spherical symmetry, the RCHB theory [39] provides a proper treatment 
of pairing correlations in the presence of the continuum and an excellent description 
for exotic nuclei [34, 28, 140, 29, 141, 38, 101, 35, 142, 143, 144, 145]. The hrst 
microscopic self-consistent description of halo in ^^Li has been achieved by using the 
RCHB theory [34]. Later, the giant halos in light and medium-heavy nuclei have been 
predicted [28, 39, 29, 144, 145]. The RCHB theory has been generalized to treat the odd 
particle system [101] and extended to including hyperons [104]. Recently, the relativistic 
Hartree-Fock-Bogoliubov (RHFB) theory in continuum has also been developed and 
used to study neutron halos [107, 114]. 

In this Section, progresses on the description of halos in spherical nuclei will be 
reviewed, including the self-consistent description of the halo in ^^Li, the giant halo and 
the extension of nuclear landscape, halo and giant halo from the non-relativistic and 
relativistic approaches, and halos from the RHFBc, etc. 

3.1. Neutron halo in 

The ground state properties of Li isotopes have been investigated by using the RCHB 
theory [34]. A satisfactory agreement with experimental values is found for the binding 
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Figure 1. Calculated and experimental density distributions in and ®Li. The 
solid line shows the result of ^^Li while the dashed line corresponds to the calculation 
for ®Li. The shaded area gives the experimental results with error bars. Taken from 
Ref. [34]. 


energies and the radii of The matter radius, showing a considerable increase 

from the nucleus ®Li to ^^Li, has been reproduced by the RCHB theory. In contrast to 
the earlier mean held calculations, e.g.. Refs. [146, 147, 96], the results given in Ref. [34] 
were obtained without any artihcial modihcations of the potential. 

In Fig. 1 the neutron density distributions of ®Li and ^^Li are shown. It is clearly 
seen that the increase of the matter radius is caused by the large neutron density 
distribution in ^^Li. Its density distribution is fully coincident with the experimental 
density shown with its error bars by the shaded area. 

By examining the mean helds for protons and neutrons and the single particle levels 
in the canonical basis, it has been shown that the halo in ^^Li is formed by Cooper- 
pairs scattered in the two levels lpi /2 and 2si/2 with the former below the Fermi surface 
and the latter in the continuum; the pairing interaction couples those levels below the 
Fermi surface with the continuum. The RCHB theory provides a much more general 
mechanism for halo: a halo appears if the Fermi level is close to the continuum and 
around the Fermi level there are single particle levels with low orbital angular momenta 
and correspondingly low centrifugal barriers. 

3.2. Giant halos and the extension of nuclear landscape 

The halo nuclei observed so far are composed of one or two halo nucleons only. The 
RCHB theory has been used to study the influence of correlations and many-body 
effects in nuclei with a larger number of neutrons distributed in the halo channel. A 
giant neutron halo has been predicted in Zr isotopes close to the neutron drip line [28]. 
It is formed by up to six neutrons scattered to weakly-bound or continuum states beyond 
the ^^^Zr core with the magic neutron number N = 82. 

Before we proceed, it is worthwhile to discuss more about the characterization 
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of a halo in medium-heavy nuclei which is one of the most relevant topics in the 
study of the halo phenomenon. The critical conditions for the appearance of a halo 
require the existence of two groups of orbitals with signihcantly different asymptotic 
slopes corresponding to the core and halo nucleons. The one-nucleon or multi-nucleon 
separation energies are close to zero, and there is a significant energy gap for the core 
and valent nucleons in the spectrum. 

Although the theoretical description of light halo nuclei is well under control, as 
discussed in Refs. [148, 149], existing definitions and tools are often too qualitative and 
the associated observables are incomplete for heavier ones. There have been lots of efforts 
in quantifying halos by examining the density profiles or the particles in the classical 
forbidden area obtained from mean field calculations [150, 151, 152]. In Ref. [148], 
a quantitative analysis method is proposed in which the halo part of the density is 
identified as the region beyond a radius tq where the core density is one order of 
magnitude smaller than the halo one. The method does not require an a priori separation 
of the density into a halo and a core part. Once rg is hxed, two criteria are introduced to 
characterize halo systems: (i) the average number of fermions participating in the halo 
and (ii) the influence of the latter on the system extension [148]. This method has been 
used to characterize halo features not only in nuclei with the spherical Hartree-Fock- 
Bogoliubov model but also in other finite many-fermion systems like atom-positron/ion- 
positronium complexes [149]. Similar attempts are highly demanded for systematic 
studies of nuclear halos within the frame of both relativistic and non-relativistic DFT. 
In this Review, we mainly analyze the halo properties based on a combined analysis of 
the following features: i) the deviation of the rms radius from the rgN^^^ law; ii) the 
large spatial extension of the density profile; and hi) the occupation of weakly bound 
and/or continuum orbitals in the canonical basis. For those weakly bound or continuum 
orbitals, we also discuss in detail their contribution to the halo. 

In the upper panel of Fig. 2, the rms radii of the protons and neutrons are presented 
for the Zr isotopes. There is a clear kink for the neutron rms radius at the magic neutron 
number N = 82. This kink can be understood by examining the single particle levels 
in the canonical basis which is shown in the lower panel of Fig. 2. Going from iV = 70 
to iV = 100, there is a big gap above the lh. 11/2 orbital. For N > 82, the neutrons 
are hlled to the levels in the continuum or weakly bound states in the order of 8 ^ 3 / 2 , 
2 / 7 / 2 , 3 pi/ 2 , 2 / 5/2 and lhg/ 2 - The neutron chemical potential is given by a dashed line. 
It approaches rapidly the continuum threshold already shortly after the magic neutron 
number N = 82 and crosses the continuum at iV = 100 for the nucleus to ^^°Zr. In 
this region the chemical potential keeps a very small but negative value for a large A 
range. This means that the additional neutrons are added with a very small, nearly 
vanishing binding energy at the edge of the continuum provided by pairing correlations. 
The total binding energies E for the isotopes above ^^^Zr are therefore almost identical. 
This has been recognized in Ref. [95]. However, the RMF calculations with the BCS 
approximation expanded on an oscillator basis is dehnitely not appropriate for exotic 
nuclei close to the continuum limit. 
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Figure 2. Upper part: Root mean square radii for neutrons and protons in Zr isotopes 
close to the neutron drip line as a function of the mass number A. Lower part: single 
particle energies for neutrons in the canonical basis. The dashed line indicates the 
chemical potential. Taken from Ref. [28]. 


From the occupation probabilities in the canonical basis, one can get the number 
of halo neutrons. It was found that there are 2 valence neutrons in ^^^Zr, 4 in ^^®Zr, 
6 in ^^®Zr, ..., roughly 5 in ^^®Zr, and roughly 6 in ^^°Zr where the neutron drip line 
is reached. Note that this does not mean that there are around six neutrons in the 
halo region of the nuclear density in ^^°Zr. We emphasize that these hgures, 2, 4, 5, or 
6, are numbers of “valence” neutrons, just as it is well recognized that there are two 
“valence” neutrons in ^^Li. With the very large neutron rms radii of these systems, one 
can estimate the number of valence neutrons which £11 in the same volume outside the 
^^^Zr core if packed with normal neutron density. This number is 24 for ^^^Zr and 34 
for ^'^‘^Zr. This phenomenon is therefore clearly a neutron halo and was called a Giant 
Halo because of the large number of particles participating in the halo formation [28]. 

It is very difficult for experimentalists to study giant halos in Zr isotopes because 
they are too heavy to be synthesized by the RIB facilities at present. Thus it is useful 
to investigate the giant halo phenomena in lighter nuclei which are easily accessible 
with available facilities. In Refs. [29, 145], ground state properties of all the even-even 
O, Ca, Ni, Zr, Sn, and Pb isotopes ranging from the proton drip line to the neutron 
drip line were investigated by using the RCHB theory with the effective interaction 
NLSH [153]. Recently, a systematic calculation for nuclei above oxygen has also been 
performed with the effective interaction PC-PKl [118], aiming to investigate the global 
impact of the continuum for the nuclear boundary; the results ranging from O to Ti 
isotopes have been published in Ref. [5]. In such a systematic study, it is found that 
the appearance of giant halos, the contribution of continuum and boundaries of nuclear 
chart are tightly connected. We note that the functional PC-PKl has proved to be very 
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successful in describing the isospin dependence of the nuclear masses [154], the Coulomb 
displacement energies between mirror nuclei [155], hssion barriers [156, 157] and nuclear 
rotations [92, 93, 91, 112], etc. 

For nuclei far from the valley of stability and with small nucleon separation energy, 
the Fermi surface is very close to the continuum threshold and the valence nucleons 
may extend over quite a wide space to form low density nuclear matter. Therefore 
the nucleon separation energies are sensitive quantities to test theory and examine the 
halo phenomena. The two-neutron separation energies for O, Ca, Ni, Zr, Sn, and Pb 
isotopes have been investigated to examine their connection with systematic behavior 
of the nuclear size [29, 144, 145]. 



N 



Figure 3. (Color online) (a) The two-neutron separation energies 5'2n calculated by 
RCHB theory [39] with NLSH [153] (Open symbols) in comparison with data available 
(solid ones) for the even O, Ca, Ni, Zr, Sn, and Pb isotopes against the neutron number 
N. (b) The root mean square (rms) neutron radii for even O, Ca, Ni, Zr, Sn, and 
Pb isotopes by RCHB calculations against the neutron number TV. The curve 
with tq = 1.139 fm has been included to guide the eye. Taken from Ref. [29] with 
modified color. 

In Fig. 3, similar to Ref. [29], the two-neutron separation energies S' 2 n and the 
neutron radii from the RCHB calculation for the even-even nuclei of the O, Ca, Ni, 
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Zr, Sn, and Pb isotope chains have been shown. The predicted curve using the 
simple empirical equation with tq = 1.139 fm normalizing to jg ^jgo 

represented in the hgure. This simple formula for agrees with the calculated neutron 
radii with exceptions in neutron rich Ca and Zr isotopes. 

The comparison between experimental and calculated 82 ^ is shown in Fig. 3. The 
experimental magic or submagic numbers N = 20, 28, and 40 are reproduced. The S' 2 n 
values for several nuclei in exotic Ca and Zr isotopes are extremely close to zero. As 
discussed in Refs. [28, 29], if taking ®°Ca or ^^^Zr respectively as a core, the valence 
neutrons will gradually occupy the loosely bound states and the continuum above the 
sub-shell of = 40 in Ca or = 82 in Zr isotopes. Such typical behavior of S' 2 n can 
be taken as an evidence of the occurrence of giant halos in Ca chain and Zr chain [28]. 

In Fig. 3, it is very interesting to see that follows the systematics well 
for stable nuclei although their proton numbers are quite different. Near the drip line, 
distinct abnormal behaviors appear at A^ = 40 in Ca isotopes and at A^ = 82 in Zr 
isotopes. As pointed in Ref. [29], nuclei having the abnormal increase correspond 
to those having small S 2 n, which provides another evidence to the emergence of giant 
halos. The increase of in exotic Ni, Sn and Pb nuclei is not as rapid as those in 
Ca and Zr chains. The regions of the abnormal increases of the neutron radii are just 
the same as those for S 2 n- Both behaviors are connected with the formation of giant 
halo [28, 29, 144, 145]. 

3.3. Halos from relativistic and non-relativistic approaches 

By using the RCHB theory, predictions of halo in the Ne [38, 158], Na [101, 158], Ca 
[29] and Zr [28] isotopes near the neutron drip line have been made. In particular, in 
Ca and Zr isotopes near the neutron drip line, the halo phenomena were addressed as 
giant halo [28, 29]. In the non-relativistic HFB approaches, Skyrme calculations with 
the parameter set SLy4 predicted the halo in Sn and Ni [152]. It was also claimed that 
the pairing gaps have an effect to reduce the halo at the neutron drip line [152, 159, 160]. 

In Refs. [161, 162], the phenomena of the giant halo and halo of the neutron-rich 
even-Ca isotopes have been investigated and compared between the frameworks of the 
RCHB and the Skyrme HFB calculations. With the two parameter sets for each of the 
RCHB and the Skyrme HFB calculations, it has been found that although the halo 
phenomena exist for Ca isotopes near the neutron drip line in both calculations, the 
halo from the Skyrme HFB calculations (SkM*) starts at a more neutron-rich nucleus 
than that from the RCHB calculations, and SLy4 does not predict a halo, because the 
drip line is much closer to the stable region than the other parameter sets. 

Figure 4 shows systematics of the rms radii of neutrons and protons fp of even- 
Ca from A^ = 30 to the two-neutron drip lines. The increase in the curvature of 
indicates the halo. For Ca, the number of nucleons in continuum, N^, of ®^“^^Ca is 0.6- 
2.2, of which the average is 1.7, in the NL-SH calculation [29]; the corresponding value 
of the SkM* calculation is always smaller than 0.5 (the average 0.27). As there are more 
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Figure 4. The neutron and proton rms radii and fp of even-Ca isotopes in RCHB 
calculations with NLSH and PKl and the Skyrme HFB calculations with SkM* and 
SLy4 in the neutron-rich region. Taken from Ref. [161]. 


than two halo neutrons, it was referred to as giant halo [28, 29]. The SkM* calculation 
predicts the halo from N = 52, and the RCHB calculations show the gradual occurrence 
of the giant halo. Apparently, the starting nucleus of the halo of SkM* corresponds to 
the extra lowering of the S 2 n- The SLy4 calculation does not predict the halo, because 
the particle-stable region ends at = 46. The large difference between fn and fp shown 
in Fig. 4 was identihed as neutron skin in the region with no halo. 

The RCHB calculations have larger neutron rms radii systematically in > 40 
than the Skyrme HFB calculations. This difference comes partially from the occupation 
of the neutron 3si/2 orbital, which causes 50% of the maximum difference in the neutron 
rms radii among the four calculations at ®®Ca. The difference in the neutron 3si/2 orbital 
is shown in Fig. 5, which illustrates the radial wave functions squared. The tail of 3si/2 
from PKl has appreciably extended distribution than the others, and as a consequence, 
the amplitude of PKl is smaller in the inner region. The number and locations of the 
nodes are the same for the four curves. 

As shown in the previous Subsection and in Ref. [28], the occupation of 3pi/2 and 
3p3/2 orbitals are responsible for the formation of giant halos in Zr isotopes. In Ref. [163], 
the Hartree-Fock-Bogoliubov (HFB) approach with Skyrme interactions SLy4 and SkI4 
has been used to study halos in Ca and Zr isotopes. It was found that the appearance 
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Figure 5. Radial wave function squared IV'rP (canonical basis) of the neutron 3si/2 
of ®®Ca. They are normalized as / dr \'4’r\‘^ = 1- Taken from Ref. [161]. 


of giant halos depends sensitively on the effective interaction adopted: SkI4 predicts a 
neutron halo in the Zr chain with A > 122 due to the weakly bound orbitals 3pi/2 and 
3p3/2- 

In the above mentioned work, the box boundary condition is adopted and the 
continuum is discretized. Although the convergence of the halo properties with respect 
to the box size has always been achieved, it is still desirable to examine the influence 
of the asymptotic behavior of the radial wave functions of those valence orbitals on the 
formation of halo or giant halo. The Skyrme HFB equations in r space have been solved 
for spherical systems with correct asymptotic boundary conditions for the continuous 
spectrum [164] or the Green’s function technique [55, 57]. In Ref. [164], it was shown that 
close to the drip line the amount of pairing correlations depends on how the continuum 
coupling is treated. In Ref. [55], it was found that, in even-even N = 86 isotones in 
the Mo-Sn region, the I = 1 broad quasiparticle resonances persist in feeling the pairing 
potential and contribute to the pairing correlation even when the widths are comparable 
with the resonance energy. In Ref. [56], the investigation of the pair correlation in the 
tail of the giant halo revealed that the asymptotic exponential behavior of the neutron 
pair condensate is dominated by nonresonant continuum states with low asymptotic 
kinetic energy. 

3.4- Halos with relativistic Hartree-Fock-Bogoliubov theory in continuum 

In spite of successes of RCHB theory, there are still a number of questions needed to 
be answered in the conventional CDFT, e.g., the contributions due to the exchange 
(Fock) terms and the pseudo-vector vr-meson. Earlier attempts to include the exchange 
terms with relativistic Hartree-Fock (RHF) method led to underbound nuclei due to the 
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missing of the meson self-interactions [165, 126]. For a long time, the RHF theory failed 
in a quantitative description of nuclear systems [166, 167, 168, 169, 170, 126, 171, 172], 
With the development of the computational facilities and numerical techniques, the 
density dependent relativistic Hartree-Fock (DDRHF) theory has been developed and 
shown signihcant improvements in a quantitative description of nuclear phenomena 
[127, 78, 113, 173, 174, 175, 176, 177] with comparable accuracy. With a number of 
adjustable parameters comparable to those of RMF Lagrangians, the DDRHF theory 
can give a equally good description of nuclear systems without dropping the Fock 
terms. Furthermore, important features like the behavior of neutron and proton effective 
masses [178] can be interpreted. 

In DDRHF theory [127], the tensor forces due to tt and p meson exchanges can 
be naturally taken into account and have brought signihcant improvements on the 
consistent description of shell evolution [173, 179, 180, 181]. The relativistic Hartree- 
Fock-Bogoliubov (RHFB) theory in continuum with density dependent meson-nucleon 
couplings is developed in Refs. [107, 114]. In the descriptions of nuclear halos, by using 
the hnite range Gogny force DIS in the pairing channel, systematic RHFBc calculations 
are performed both for stable and weakly bound nuclei. It is demonstrated that an 
appropriate description of both mean held and pairing ehects can be obtained within 
RHFBc theory [107, 114]. 

In Fig. 6, the single-neutron separation energies Sn of Sn isotopes from ^°^Sn to 
^^®Sn (left panels) and the single-proton separation energies Sp oi N = 82 isotones from 
^^°Cd to ^^^Lu (right panels) are given. The odd-even diherences on the single-nucleon 
separation energies rehect the ehects of the pairing correlations. In Fig. 6, RHFBc with 
PKAl [113] and PKOl [127] as well as RHB with DD-ME2 [182] present comparable and 
satished quantitative agreements with the data for both isotopes and isotones. From 
Fig. 6, one can hnd some systematics in the results. On the neutron rich side, i.e., 
after ^^^Sn for Sn isotopes and before ^ 64 Gfd for N = 82 isotones, PKAl shows better 
agreements than PKOl and DD-ME2. This is mainly due to the improvements in the 
shell structures with exchange terms included, especially the shell closures beyond 50 
and 82, i.e., 58 and 92, which may change the level density and then inhuence the pairing 
ehects. On the proton rich side, these three approaches provide comparable agreement 
with the data. 

The Ce isotopes have been chosen to demonstrate the nuclear halo phenomenon 
within the RHFBc theory [107, 114]. In Fig. 7 the nuclear matter distributions (left 
panels) and neutron canonical single particle conhgurations (right panels) for Ce isotopes 
close to the drip line are shown. The neutron drip line calculated here with PKAl [113] 
is A^ = 140 whereas the calculations with PKOl [127] and DD-ME2 [182] predict a 
shorter one as N = 126. As shown in Fig. 7(a), the neutron densities become more and 
more diffuse after the isotope ^®®Ce (A^ = 128), a direct and distinct evidence of halo 
occurrence. 

From Fig. 7(b) one can see that such extremely extensive matter distribution, e.g. 
in ^®®Ce, is mainly due to the low-/ states, namely the halo orbitals i^4si/2, i^ 8 d 5 j 2 and 
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Figure 6. (Color online) Single-nucleon separation energies along Sn isotopic {Sn- left 
panels) and N = 82 isotonic {Sp: right panels) chains. The results are calculated by 
RHFB with PKAl [113], PKOl [127], and by RHB with DD-ME2 [182], in comparison 
to the experimental data [183]. Taken from Ref. [114]. 


u3ds/2- From the occupations of the halo orbitals iVhaio in Fig. 7(d), the halos were 
predicted in ^®®Ce, ^®^Ce and ^®°Ce and giant halos may appear in ^®^Ce, ^®^Ce, ^®®Ce 
and ^®®Ce because more than two neutrons are occupying the halo orbitals. Similar 
conclusions can be obtained from the neutron numbers lying beyond the sphere with 
the radius r = 10 fm [A7->io fm in Fig. 7(d)], which is large enough (the neutron matter 
radius = 6.2 fm in ^®®Ce) for halos. Even extending to r = 16 fm, which is sometime 
taken as the radial cut-off in the calculations of stable nuclei, there are still some amount 
of neutrons lying beyond this sphere in the isotopes from ^®^Ce to ^®®Ce. 

In Fig. 7(c), one Ends that the halo orbitals (i^4si/2, i'^d ^/2 and i'3d3/2) are located 
around the particle continuum threshold where they are gradually occupied. For the 
isotopes beyond ^^'^Ce {N = 126), the Fermi levels (in open circles) approach the 
continuum threshold rather closely such that the stability of these halo isotopes becomes 
sensitive to pairing effects. Nearby the low-/ states, there are the high-/ states v2g^i2 
and u2gjj2. Because of the relatively large centrifugal barrier for ^f-orbitals they do not 
contribute much to the diffuse neutron distributions. Nevertheless, the existence of the 
high-/ states nearby halo orbitals is still particulary significant, because it leads to a 
rather high level density around the Fermi surface, and evidently the pairing effects are 
enhanced to stabilize the halo isotopes. 

Evidence for the existence of a halo can also be investigated by the systematic 
behavior of nuclear bulk properties such as radii. In Fig. 8, the isospin dependence of 
the neutron skin thickness (r„ — Vp) calculated in density dependent RHFB theory in 
continuum using PKAl for Ca, Ni, Zr, Sn, and Ce isotopes are given [107]. Taking the 
stable nuclei as references (shown as the dashed lines), continuously growing deviations 
are found in the Ca, Zr, and Ce chains until the neutron drip line. This can be considered 
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Figure 7. (Color online) (a) Neutron and proton densities, (b) the relative 
contributions of different orbitals to the full neutron density in ^®®Ce, (c) neutron 
canonical single particle energies occupation probability (in x-error bars) and Fermi 
energy Ep (in open circles), and (d) neutron numbers filling in the halo orbitals 
4 si/ 2, 3(i5/2, and 3^3/2 (-Nhaio), and the ones lying beyond the spheres with the 
radii r = 10, 11, 12, 13, 14, 15, 16 fm, respectively Nr>w fm, Nr>ii fm, Nr>i 2 fm, 
Nr>i 3 fm, fm, fm and fm- The results are calculated by RHFB 

with PKAl [113] plus the Gogny pairing force DIS [124]. The spherical box radius is 
adopted as i? = 28 fm. Taken from Ref. [107]. 


as evidence of a halo. Despite the deviations in the mid-region, both Ni and Sn nuclei 
show a similar isospin dependence in both stable and exotic nuclei regions, which indicate 
a neutron skin since the growth of a halo is interrupted. Compared with Ni, the Sn 
isotopes show even a much weaker skin effect. 

The systematics for the Ca and Zr chains from the RHFBc calculations are similar 
to those from the RCHB calculations [29]. For the Ce isotopes, RHFBc calculations 
with PKAl show clear evidence for the existence of halo structures. The calculations of 
RHFBc with PKOl [127] and RCHB with DD-ME2 [182] predict that the isotopic chain 
ends at = 126, before the halo occurrence as predicted by RHFB with PKAl. This 
deviation is interpreted by the shell structure evolution in Fig. 7(c), where much weaker 
shell are predicted with PKAl in comparison with PKOl and DD-ME2. As shown in 
Fig. 7(c), the neutron shell gap (A^ = 126) between Z/IA 3/2 and states is close to 

the particle continuum threshold, which might essentially influence the stability of the 
drip line isotopes. 

From Ref. [107], one may conclude that the Fock term itself has little direct impact 
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N-Z 


Figure 8. (Color online) Neutron skin thickness r„ — Vp for Ca, Zr, Ni, Sn, and Ce 
as a function of the isospin (TV — Z). The r„ and Vp are respectively the neutron and 
proton root mean square radii. The results are calculated by RHFB with PKAl [113] 
plus the paring force DIS [124], Taken from Ref. [107]. 


on the formation of a nuclear halo. However, the Fock term may influence the shell 
structure, particularly those orbitals close to the Fermi surface in drip-line nuclei, thus 
affecting indirectly the formation of a potential halo. 

3.5. Halos in exotic hyper nuclei 

Motivated by knowledge of A-N interaction and understanding on giant halo [28, 29, 145], 
it is interesting to study the possible appearances of halos in exotic hyper nuclei. 

Two neutron separation energies S 2 n for normal nuclei, single-A hyper nuclei and 
double-A hyper nuclei of Ca isotopes, labelled by A = 0, A = 1 and A = 2, respectively, 
from the proton drip line to neutron drip line were investigated [104]. It was found that 
one or two A hyperons lower the Fermi level due to the attractive A-N potential but 
keep the neutron shell structure unchanged. Therefore, the neutron drip line is pushed 
outside from TV = 52 in Ca isotope chain to iV = 54 in hyper isotope chain. Meanwhile, 
giant halos due to pairing correlation and the contribution from the continuum still exist 
in Ca hyper nuclei similar to that in Ca isotopes [29]. This is a slight but rewarding 
step for exploring the limit of drip line nuclei on the basis of the giant halo. For more 
details, see Ref. [104]. 

Apart from neutron halos in hyper nuclei, as A hyperon is less bound than the 
corresponding nucleon in nuclei, it is worth investigating the existence of hyperon halos. 

In Ref. [184, 185], a A halo was tentatively suggested in a RMF plus BCS model. 
However, since the normal BCS method will have unphysical solution by involving 
baryon gas [39], such prediction needs further study. In Ref. [103], hyper carbon isotopes 
are studied by the RCHB theory. By investigating the baryon density distributions in 
hyper nuclei j^^C, 2 aC, and g^C, it was found that with up to two A hyperons added 
to the core ^^C, the nucleon density distributions remain the same and hyperon density 
distributions at the tail are comparable with those of the nucleons. An intriguing 
phenomenon appears in g^C where the hyperon density distribution has a long tail 








CONTENTS 


27 


extended far outside of its core 2 aC. This is a signature of hyperon halo due to the 
weakly bound state IP 3/2 in sa^ which has a small hyperon separation energy and a 
density distribution with long tail [103]. 


4. Surface diffuseness 

From the mean held point of view, the properties of nucleons in an atomic nucleus 
are determined by the mean potential which stems from their interaction with the other 
nucleons and is determined by nucleon densities in the self-consistent models. Therefore, 
the study of the isospin dependence of the mean potential and the density distribution, 
both become highly diffuse near the particle drip line, is crucial to understand exotic 
nuclei. For hnite nuclei, the diffuseness of the nuclear surface provides a measure of the 
thickness of the surface region. The surface diffuseness is intimately related to the spin- 
orbit splitting [142], the evolution of single particle shell structure [28], the pseudospin 
symmetry [76, 77, 186] and spin symmetry [83, 186], and the nuclear surface energy and 
thus nuclear masses [187]. In this Section, we will review recent progresses on nuclear 
surface diffuseness and the corresponding impact on nuclear shell structure and size. 


4-1. Spin-orbit splitting and shell structure 


The RCHB theory has been used to systematically study the isospin dependence of the 
mean potential for the Sn isotopes from proton drip line to neutron drip line and the 
spin-orbit splitting for the whole isotope chain was examined in details [142]. As an 
example, the spin-orbit splitting 


Els — 


Eij=i_i/2 — F'/j=Z+l/2 
2 / + 1 


(60) 


versus the binding energy: 

^ (/ + + lEij=i+i/2 

“ 2/ +1 ^ ^ 
in "^Sn with A = 110,120, • • •, 170 are given in Fig. 9 for the neutron spin-orbit 
partners ( 1 ^ 3 / 2 , 1 ^ 5 / 2 ), ( 15 - 7 / 2 , l 5 ' 9 / 2 ), ( 1 Ai/ 2 , IA 3 / 2 ), (lpi/ 2 , IP 3 / 2 ), ( 1 / 5 / 2 , 1 / 772 ) and 
( 1 / 19 / 2 , 1 / 111 / 2 ), and the proton spin-orbit partners ( 1 ^ 3 / 2 , 1 ^ 5 / 2 ) and (I/ 5 / 2 , 1 / 7 / 2 )- 

It is very interesting to see that the spin-orbit splittings for the neutron and proton 
are very close to each other, at least for (lc/3/2, 1 ^ 5 / 2 ) and (I/ 5 / 2 , 1 / 7 / 2 ) cases. The 
splitting decreases monotonically from the proton drip line to the neutron drip line. 
To understand its mechanism, it is very helpful to examine the origin of the spin-orbit 
splitting in the Dirac equation. 

For the Dirac nucleon moving in scalar and vector potentials, its equation of motion 
could be de-coup led and reduced for the upper component and the lower component, 
respectively. For the spin-orbit splitting, the Dirac equation can be reduced for the 
upper component as following: 

|. 1 d{2M — Tjo +'Es) 

^d^ E + 2M -Eq + Es d/ d^^ * 
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Figure 9 . The neutron spin-orbit splitting Eis = — ^'ij=z-i-i/2)/(2^ + 1 ) 

versus the mass number A in Sn isotopes for neutron (1013/2,1^5/2)1 (l<?7/2, lff9/2), 
(1*11/2,1*13/2), (lpi/2, IP3/2), (1/5/2,1/7/2) and (1/19/2,11*11/2) orbitals and proton 
(l(i3/2, l(i5/2) and (1/5/211/7/2) orbitals, respectively. Taken from Ref. [ 142 ]. 


k(1 A- k) 


K d{2M — So + S 5 ) 




E + 2M — Tq + Tis r dr 

= - (E + 2M - So + J:s){E - So - ^s)G^Hr)- 
The spin-orbit splitting is dne to the corresponding spin-orbit potential 
1 K d(2M — Sq -|- S5) 


(62) 


(63) 


E -b 2M — So -b S 5 r dr 

with some proper normalization factor. The spin-orbit splitting is energy dependent 

and depends on the derivative of the potential 2M — Sq -b S5. Therefore the so-called 

K d{2M — So -b S5) 


spin-orbit potential, Vis = 


r dr 


is introdnced. 


The derivatives of the nentron potentials So(r) — S5(r) for ^^°Sn, and ^'^°Sn, 

are given in the npper panel of Fig. 10. For both proton and nentron, the derivatives 
of the potentials So(r) — S5(r) are almost the same. As the potential So — S5 is a 
big qnantity (~ 700 MeV), the isospin dependence in the spin-orbit potential could be 
neglected. That is the reason why the spin-orbit splitting for the neutron and proton is 
very close to each other in Fig. 9. From ^^°Sn to the amplitude of the derivative 

for So — S 5 decreases monotonically due to the surface diffuseness. 

The decline of the spin-orbit splitting is from the diffuseness of the potential or 
the outwards tendency of the potential. The neutron and proton potentials So -b S5 
are given in Fig. 11 for "^Sn with A = 110,120, • • •, 170. It is seen that the depth of 
the neutron potential decreases monotonically from the proton drip line to the neutron 
drip line and the surface moves outwards. The inserted hgures give the radii Rq at 
which So -b S5 = —10 MeV as functions of the mass number. With the increase of 
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Figure 10 . Upper: The derivative of the neutron potentials Eo(?') ~ 5]s(r) for "^Sn 
with A = 110,120, • • •, 170; Lower: The neutron potentials So(r) — Ss(r) for "^Sn with 
A = 110,120, • • •, 170. In order to examine the surface diffuseness more clearly, the 
radii Rq at which Eo(i?o) — ^s{Ro) = 100 MeV has been given as an inserted figure. 
Taken from Ref. [142]. 


the neutron number, the depth of the proton potential increases monotonically but is 
pushed towards outside as well due to the proton-neutron interaction. 

Similar neutron potentials Sq — for "^Sn with A = 110,120, • • •, 170 are given in 
the lower panel of Fig. 10 and its inserted hgure gives the radii Rq at which Sq —F 5 = 100 
MeV as a function of the mass number. As seen above, the spin-orbit splitting is related 
with the derivative of the potential Sq — S 5 . The surface diffuseness happens for both 
the vector and scalar potential as shown in both Sq — S 5 and Sq -|- 

Similar RCHB calculations have been carried out for all the nuclei in Na isotopes 
with mass number A ranging from 17 to 45 with effective interaction NLSH [28]. Apart 
from the diffuseness of the potential, the shell evolution with neutron number has been 
examined there as well. In Fig. 12, the microscopic structure of the single particle 
energies in the canonical basis [188, 39] is given. In the left panel of Fig. 12, the single 
particle levels in the canonical basis for the Na isotopes with an even neutron number are 
shown. The level density becomes dense with the neutron number due to the decrease 
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Figure 11. Left: The neutron potentials So(''’) + '^sO for ^Sn with A = 
110,120, • • •, 170. In order to examine the surface diffuseness more clearly, the radii 
i?o at which So(i?o) + ^s{Ro) = —10 MeV has been given as an inserted figure; 
Right: The proton potentials So(?') + Ss(r) for ^Sn with A = 110,120, • • •, 170. 
In order to examine the surface diffuseness more clearly, the radii Rq at which 
So(i?o) + Ss(i?o) = —10 MeV has been given as an inserted hgure. Taken from 
Ref. [142]. 

of the spin-orbit potential or diffuseness of the potential. Going from A = 19 to A = 45, 
it can be observed a big gap above the N = 8, N = 20 major shell, and iV = 14 sub¬ 
shell. The N = 28 shell for stable nuclei fails to appear, as the 2p3/2 and 2pi/2 come so 
close to 1 / 7 / 2 - When N > 20, the neutrons are hlled to the levels in the continuum or 
weakly bound states in the order of I/7/2, 2 ^ 3 / 2 , 2 pi/ 2 , and 1 / 5 / 2 - In the right part, the 
occupation probabilities in the canonical basis of all the neutron levels below E = 10 
MeV have been given for ^^Na to show how the levels are hlled in nuclei near the drip 
line. The importance of careful treatment of the pairing correlation and the scattering 
of particle pairs to higher lying levels are noted in the hgure. 

4 . 2 . Interaction cross section 

The neutron and proton distributions and radii are fundamental quantities to dehne the 
halo phenomena. However, the radius for unstable nucleus is normally extracted from 
the interaction cross sections by assuming an empirical form (e.g., a Gaussian function) 
for the density distribution [24]. It would be better to use the density distribution 
from the self-consistent models and calculate the interaction cross sections directly from 
Glauber model. 

Such systematic calculations have been carried out for all the nuclei in Na isotopes 
with mass number A ranging from 17 to 45 with the RGHB theory and the ehective 
interaction NLSH [101]. The calculated binding energies E and the interaction cross 
sections with the Glauber model are shown in Fig. 13. The calculated binding energies 
E are compared with the empirical values [189]. The proton (neutron) drip line nucleus 
has been predicted to be ^®Na (^^Na). The difference between the calculations and 
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A ==Na, 

Figure 12. Left part: Single particle energies for neutrons in the canonical basis as a 
function of the mass number for Na isotopes. The dashed line indicates the chemical 
potential. Right part: The occupation probabilities in the canonical basis for ^^Na. 
Taken from Ref. [101]. 


the empirical values for the stable isotopes is from the deformation, which has been 
neglected in Ref. [101]. 

As the interaction cross sections ai is well reproduced, the density distributions of 
the whole isotopes from the RCHB theory have been examined and the relation between 
the development of halo and shell effect has been studied in detail [101]. It was found 
that the neutron N = 28 shell closure fades due to the lowering of the 2p orbitals. 
The tail in the nuclear density prohle depends strongly on the shell structure and the 
development of a halo was connected with the changes in the occupation in the next 
shell or sub-shell close to the continuum limit [101]. 

Similar calculations have been performed for C, N, O and F isotopes up to the 
neutron drip line by the self-consistent RCHB theory [143]. The agreement between the 
calculated one neutron separation energies Sn and the available experimental values are 
similar as Na isotopes. A Glauber model calculation for the total charge-changing cross 
section has been carried out with the density distribution obtained from the RCHB 
theory. The measured cross sections for studied nuclei with ^^C as a target are nicely 
reproduced. An important conclusion was found that, contrary to the usual impression, 
the proton density distribution is less sensitive to the neutron number along the isotope 
chain. Instead it is almost unchanged from stability to the neutron drip line. Such 
calculations are quite useful in extracting both the proton and neutron distributions 
inside the nucleus and in dehning the neutron/proton skin or halo. 
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16 20 24 28 32 36 40 44 

A 

Figure 13 . Upper part: The interaction cross sections ai of "^Na isotopes on a 
carbon target at 950A MeV: the open circles are the result of RCHB calculation and 
the available experimental data {A =20-23, 25-32) are given by solid circles with their 
error bars. The dashed line is a simple extrapolation based on the RCHB calculation 
for Lower part: Binding energies for Na isotopes, the convention is the same 

as the upper part, but the RCHB result for particle unstable isotopes are indicated by 
triangle. Taken from Ref. [101]. 


4-3. Surface diffuseness correction in global mass formula 

For finite nuclei, the diffuseness of nuclear surface is an important degree-of-freedom 
in the calculations of nuclear masses [187]. The density distributions of most stable 
nuclei are of the “neutron skin-type”, with a typical value around 0.5 fm for the surface 
diffuseness. For nuclei near the neutron drip line, halos and giant halos may develop and 
the neutrons may distribute more extensive spatially than the protons, which implies the 
enhanced neutron surface diffuseness for these extremely neutron-rich nuclei. However, 
in global mass calculations, the surface diffuseness of exotic nuclei near the drip lines 
has not been properly considered. Recently, the surface diffuseness from the RCHB 
theory and its impact on the symmetry energy and shell correction have been studied in 
Ref. [187] and the accuracy of the global mass formula has been considerably improved. 

Inspired by the Skyrme energy-density functional, a macroscopic-microscopic mass 
formula, the Weizsacker-Skyrme (WS) formula [190, 191, 192], was proposed with a rms 
deviation around 336 keV with respect to the 2149 measured masses [183] in 2003 Atomic 
Mass Evaluation (AME). The WS models provide the best accuracy for nuclear masses 
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Figure 14 . (Color online) (a) Difference between measured and the WS4 calculated 
masses, (b) Difference between WS3 and WS4 calculated masses. The squares and 
crosses denote the nuclei first appearing in AME2012 and the nuclei with |/ —/o| > 0.1, 
respectively. The smooth and the zigzag curves denote the /3-stability line from the 
Green’s formula and the neutron drip line from the WS4 formula, respectively. Taken 
from Ref. [187]. 


in several mass regions [193, 194] and has been used in r-process simulations [195, 196]. 

In the WS formula, an axially deformed Woods-Saxon potential, with a contant 
surface diffuseness parameter a for all nuclei, was used to obtain the single particle 
levels of nuclei. In Ref. [187], by taking into account the surface diffuseness correction 
for unstable nuclei, a new global mass formula, WS4, has been proposed. The surface 
diffuseness a of the Woods-Saxon potential was given by a = ao{l + 2e6q), with 
e = {I — IqY — the correction factor for surface diffuseness, oq the diffuseness of the 
Woods-Saxon potential, Jq = 0.4y4/(R-|-200) the isospin asymmetry of the nuclei along 
the /9-stability line described by Green’s formula, and 6q = 1 for neutrons (protons) in 
the nuclei with / > Jq (/ < Jq) and Sq = 0 for other cases. This means that the surface 
diffuseness of neutron distribution is larger than that of protons at the neutron-rich 
side and smaller than that of protons at the proton-rich side. The rms deviation with 
the available mass data [197] falls to 298 keV, crossing for the hrst time the 0.3 MeV 
accuracy threshold for mass formulas (models) within the mean field framework. 

In Fig. 14(a), it was shown the deviations of the calculated masses from the 
experimental values. For all the 2353 nuclei with Z and N > 8, the deviation is within 
1.23 MeV. In Fig. 14(b), the difference between the WS3 (constant surface diffuseness Oq) 
and WS4 calculated masses. For most nuclei, the results of WS3 and WS4 are consistent 
in general (with deviations smaller than one MeV). For nuclei near the neutron drip line, 
the masses given by WS4 are larger than the results of WS3 by several MeV. This is due 
to the enhancement of the nuclear symmetry energy coming from the surface diffuseness 
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effect in the extremely neutron-rich nuclei. 

It is worthwhile to note that with an accuracy of 258 keV for all the available 
neutron separation energies and of 237 keV for the alpha-decay Q-values of superheavy 
nuclei, the proposed mass formula will be important not only for the reliable description 
of the r-process of nucleosynthesis but also for the study of the synthesis of superheavy 
nuclei. 


5. Pairing correlations and nnclear size 


In Ref. [159], it is shown that the neutron ground state Hartree-Fock (HF) asymptotic 
density equals 

p(r) oc exp (—2pr)/r^, (64) 


with /i = y/—2mek/h and Sk the single particle energy of the least bound I = 0 neutron. 
The corresponding mean square radius deduced from the asymptotic solution is 


(r^)HF oc 


(65) 


2 m \ ek \’ 

which diverges in the limit Sk —)■ 0. 

In the presence of pairing, the mean square radius deduced from the asymptotic 
Hartree-Fock-Bogoliubov (HFB) density is 


{r )hfb oc 


2m{Ek — A) 


( 66 ) 


with the lowest discrete quasiparticle energy Ek = -s/(sk — A)^ -|- A^. If the pairing 
gap Ak is hnite, the radius will never diverge in the limit of small separation energy 
Sk — A —y 0. 

The asymptotic HF and HFB densities characterized by / = 0 orbitals were 
compared in Ref. [159] and it has been emphasized that pairing correlations reduce 
the nuclear size and an extreme halo with inhnite radius cannot be formed in superfluid 
nuclear systems. Here pairing correlations act against the formation of an inhnite radius. 
Therefore this effect was called “pairing anti-halo effect” in Ref. [159]. 

In fact, the mechanism leading to a halo in Eq. (65), i.e., the valence nucleons 
occupy the weakly-bound orbital below and very close to the continuum threshold, has 
been used in early interpretation for nuclear halo phenomena [146, 147, 96]. The limiting 
condition —)■ 0 and the radius deduced from this I = 0 orbital alone correspond to an 
extremely ideal situation, which is difficult to be found in real nuclei. 

In Ref. [198], the inhuences of pairing correlations on the nuclear size and on the 
development of a nuclear halo were studied in details by the self-consistent RCHB 
theory [39]. In order to simplify the problem, neutron-rich nuclei with the neutron 
Fermi surface below, between and above two weakly-bound 2p levels are investigated 
with a hxed Wood-Saxon potential, which is followed by the self-consistent calculation 
on the well studied neutron halos in ^^Li and ^^Ne with low-/ orbitals in the continuum. 
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5.1. Fermi surface and nuclear size 

In order to investigate the impact of nnclear binding on nnclear size for varying 
pairing potential, a spherical Woods-Saxon shape mean held potential htted to the 
corresponding neutron potentials resulting from a self-consistent calculation for '^^Mg is 
adopted [198]. 

The total neutron radius of the nucleus is given by 

^ (67) 

nlj 

which is determined by the rms radii of the orbitals {r)nij and the corresponding 
occupation probabilities 



A (MeV) A (MeV) A (MeV) A (MeV) 


Figure 15 . (Color online) Single neutron levels 2p and If near the Fermi surface in 
the canonical basis, the rms radii {r)nij, occupation probabilities contributions 
X {r)nij to the neutron rms radii, and the total neutron rms radii i?jv for 
N = 28, 30, 32, and 34 as a function of average pairing gap A. The Fermi surfaces 
for the neutron numbers are plotted as dashed lines. Taken from Ref. [198]. 

In Fig. 15, the single neutron levels Ek in the canonical basis [188], the rms radii 
{'>^)nij, the occupation probabilities ^nlj^ the contributions v'^ij x {r)nij to the neutron 
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rms radius for the orbitals near the Fermi surface {2p and 1/), and the total neutron 
rms radii Rjq for four neutron-rich Mg isotopes with N = 28, 30, 32, and 34 are plotted 
as functions of the pairing gap 

+ ( 68 ) 

nlj 

Here are the diagonal matrix elements of the pairing field in the canonical basis. 

The single neutron levels in canonical basis are obtained by solving the 
corresponding RCHB equations for various pairing strength. In principle, the energy 
levels Enij = hnn depend on A. However, since this solution was obtained for fixed 
WS potentials, as shown in the upper panels of Fig. 15, these states remain almost 
unchanged with increasing pairing correlations. 

For Mg isotopes with N = 28, 30, 32, and 34, the neutron Fermi surface (shown as 
dashed line in upper panels of Fig. 15) is located just below the two weakly-bound 2p 
levels, between them, crossing the 2 pi/2 level, and above the 2 p levels. As the pairing 
strength increases, more neutrons are scattered from occupied levels in the Fermi sea to 
empty levels above the Fermi level, and therefore the Fermi surface is raising. 

In Fig. 15, the rms radii of the 2p orbitals are much larger than those of the 
If orbitals due to the lower centrifugal barrier. With increasing pairing correlations, 
the rms radius decreases for 2p and 1/ levels. This is the so-called “pairing anti-halo 
effect” discussed in Ref. [159]. However, this effect concerns only the radii {r)nij of the 
individual orbitals. The total neutron radius of the nucleus in Eq. (67) is determined 
by the rms radii of the orbitals and the occupation probabilities, which depend strongly 
on the pairing correlations. 

As shown in Fig. 15, for N = 28 the neutron Fermi surface is just below the two 
weakly-bound 2p levels. Without pairing correlations, the occupation probability is 1.0 
for the 1 / 7/2 orbital, and it vanishes for the 2 p and the I/5/2 orbitals. As the pairing 
strength increases, the neutrons on I/7/2 orbital are scattered to 2 p and I/5/2 orbitals 
which have much larger rms radii. Therefore, the contributions to the total neutron 
rms radius from 2 p and I/5/2 states grow more than the contribution from I/7/2 state 
decreases. As a result, the total neutron rms radius increases monotonically, where the 
contributions from 2 p orbitals play a dominant role. 

For A = 30 in Fig. 15 , the 2pi/2 and the I/5/2 orbitals are empty for zero pairing, 
while the I/7/2 orbital is fully occupied and the 2 ^ 3/2 orbital is half occupied. With 
increasing pairing strength, the neutrons in the I/7/2 orbital are scattered to the 2 p 
orbitals which contribute strongly to the neutron rms radius. Therefore, the total 
neutron rms radius Rn increases with the growing pairing strength up to A = 0.9 
MeV. Increasing the pairing strength further, the neutrons are continuously to be 
scattered from the I/7/2 to the 2 pi /2 state, while the neutron number in the 2 p ^/2 orbital 
decreases slightly and the resonance state I/5/2 begins to be occupied. Together with 
the monotonic decreasing rms radii of the individual orbitals, the increase resulting from 
the 2 pi /2 and I/5/2 states is smaller than the decreases of the contributions from the 
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2^3/2 and I/7/2 states. Therefore, for very strong pairing correlations, the total neutron 
rms radius dually decreases. As it is clearly seen this decrease is not caused by the 
pairing anti-halo effect, because the radius of the I/7/2 state stays completely constant. 
It has its origin in the reoccupation caused by pairing. 

For = 32 in Fig. 15, the I/7/2 and 2^3/2 orbitals are fully occupied without 
pairing correlations, while the 2pi/2 and I/5/2 orbitals are empty. As the pairing strength 
increases, the neutrons in the I/7/2 and 2^3/2 orbitals are scattered to the 2pi/2 orbital, 
whose rms radius is larger than the ones of I/7/2 and 2^3/2 levels. The contribution 
of the 2pi/2 orbitals causes an increase in the total neutron rms radius up to a pairing 
gap A = 0.9 MeV. As the pairing strength continues to grow beyond this value, the 
neutrons in the 2^3/2 state begin to be scattered to the I/5/2 state which provides 
a smaller contribution than the 2^3/2 state, while the occupation probability remains 
almost constant for the lpi/2 state. Together with the decreasing rms radius of the 
individual orbitals, the total neutron rms radius dually decreases for very large pairing. 

For N = 34, the two weakly-bound 2p levels are fully occupied for zero pairing, 
and the Fermi surface is just above the weakly-bound 2p levels. With increasing pairing 
strength, the neutrons in the I/7/2 and 2p orbitals are scattered to the I/5/2 orbital 
with a much smaller rms radius than the 2p orbitals. The contributions to the total 
neutron rms radius from the I/7/2 and 2p orbitals decrease more than the increase of 
the contribution from I/5/2 orbital, and therefore the total neutron rms radius dually 
decreases monotonically. Here the two edects, decreasing rms radius of the individual 
orbitals and changes of the occupation probabilities, act in the same direction. The 
decrease of the total neutron radius with increasing pairing correlations is not only 
from the decreasing rms radii of the individual orbitals but also from the change of the 
occupation probabilities. 

It is clear that pairing correlations can change the rms radii of the individual 
weakly-bound orbitals and simultaneously the corresponding occupation probabilities. 
As a result, they have an induence on the nuclear rms radius and the total nuclear 
size. Furthermore, the contributions of the weakly-bound 2p orbitals to the nuclear rms 
radius play an important role. 

5.2. Crucial low I orbitals 

Following the investigation on the impact of nuclear binding on nuclear size, the next 
question is the impact of orbital binding on the nuclear size as a function of the pairing 
potential. In Ref. [198], this problem has been investigated for the orbitals 2p as well- 
bound, weakly-bound, around the continuum threshold, and in the continuum cases. 

In Fig. 16(a), the single neutron levels 1^3/2, I/7/2, 2^3/2, 2.pi/2, and I/5/2 for 
A = 28 in WS potentials with varying depth such that £2^1/2 = —5.0, — 0.5, 0.0, and 
-1-0.5 MeV are shown. In Fig. 16(b), for each £2^1/2; ^^e change of the total neutron rms 
radius with respect to that of zero pairing, SRn = Rat (A) — Rat (A = 0), is plotted as 
a function of the pairing gap. 


CONTENTS 


38 



S2PI/2 (MeV) 



A (MeV) 


Figure 16 . (a) Single neutron levels 1^3/2, I/7/2, 2^3/2, 2pi/2, and I/ 5/2 in WS 

potentials of varying depth with e 2 pi /2 = —5.0, — 0.5, 0.0, and 0.5 MeV. (b) Changes 
of the total neutron rms radii 5Rn = Rn{N) — Rn{^ = 0 ) for the various cases of 
£^ 2 pi /2 as a function of average pairing gap A. The shadow area represents the results 
calculated with different box sizes in the WS potential for £ 2 ^ 1/2 = 0.5 MeV. Taken 
from Ref. [198]. 


For e 2 pi /2 = —5.0 MeV which corresponds to a well bound nucleus, the neutron rms 
radius Ri^{N) increases slightly with the pairing gap increasing from 0.0 to 1.5 MeV. For 
^ 2 pi /2 = —0.5, 0.0, + 0.5 MeV which correspond to a weakly-bound nucleus, increases 
by more than 0.1 fm are observed. It has to be noticed that for £ 2 pi /2 = +0.5 MeV, 
both levels 2pi/2 and 2 ^ 3/2 are in the continuum. If the pairing correlation is switched 
on, the system is basically not bound. The eigenvalue problem is solved in a hnite box 
and the solution depends on the box size. As an illustration, when the box size changes 
from 20 to 25 fm, the neutron rms radius calculated with a box radius of 25 fm is about 
0.05 fm larger than that obtained with with a radius of 20 fm box for A = 1.5 MeV as 
shown as shadowed region in Fig. 16(b). 

For the discussion on the influence of pairing correlations on the nuclear size so 
far, the mean helds are hxed in the form of spherical WS potentials, and only the 
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pairing strength and the Fermi level for N = 28, 30, 32, 34 are considered. Here, 
with the changing pairing strength and the Fermi level, the single neutron levels stay 
almost constant except slight modihcation by the pairing held. In Ref. [198], fully self- 
consistent RCHB calculation for neutron-rich nuclei '^o,42,44,46]y[g been performed, it 
should be emphasized that qualitatively the effect of pairing correlations on the nuclear 
size agrees with that found for hxed WS potentials. 
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Figure 17 . (Color online) Same as Fig. 15, but for and ^^Ne. Taken from 
Ref. [198]. 
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5.3. Pairing correlations and halos 

In order to investigate the effect of pairing correlations on the development of a nuclear 
halo, self-consistent calculations are performed with different pairing strength for the 
well known neutron halo nuclei ^^Li [34] and ^^Ne [38], where the low-/ orbitals are in the 
continuum [198]. Single neutron levels in the canonical basis, the rms radii {r)nij, the 
occupation probabilities the contributions to the neutron rms radius v'^ij x {r)nij, 
and the total neutron rms radii are plotted in Fig. 17 as a function of the average pairing 
gap for the nuclei ^^Li (panel a) and ^^Ne (panel b). 
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For in Fig. 17(a), a halo is developed by scattering the neutrons to the 
2 si /2 orbitals in continuum and the increasing pairing correlations will promote the 
development of this halo. 

For ^^Ne in Fig. 17(b), the 2p and I/7/2 orbitals are in the continuum. Without 
pairing correlations, the last two neutrons occupy the 2 ^ 3/2 orbital and therefore this 
nucleus is unbound and the rms radius depends on the box size. Increasing the pairing 
strength, the Fermi surface comes down and becomes negative when the pairing gap 
gets larger than 1.0 MeV. It shows that for the halo nuclei such as ^^Ne, where the low-/ 
orbitals are in the continuum and unoccupied without pairing correlation, increasing 
pairing correlations will decrease the size of the nuclear halo, but nonetheless pairing 
plays an essential role: Without pairing the nucleus would stay unbound. 

It is concluded that the weakly-bounded orbitals with low orbital angular momenta 
/ in the neighborhood of the threshold to the continuum play an important role [198]. 
The pairing correlations have a two-fold influence on the total nuclear rms radius and 
the nuclear size. First, they can change the rms radii of individual weakly bound 
orbitals and, second, they can change the occupation probabilities of these orbitals. 
With increasing pairing strength, the individual rms radii are reduced (pairing anti¬ 
halo effect). Meanwhile pairing changes also the occupation probabilities. The total 
nuclear size is determined by the competition between this two effects. 

In Ref. [198], the strength of the pairing correlations was an external variable 
parameter. In realistic nuclei, the strength of pairing depends on the level density 
in the vicinity of the Fermi surface. Therefore, even orbitals with high /-values and 
large centrifugal barriers can contribute indirectly to changes of the nuclear size [38]. If 
they come close to the Fermi surface, they enhance pairing correlations and influence at 
the same time the occupation of low-/ orbitals close to the continuum limit and therefore 
the nuclear size. 

6. Deformed halo and shape decoupling 

Most open shell nuclei are deformed. The interplay between deformation and weak 
binding raises interesting questions, such as whether or not there exist halos in deformed 
nuclei and, if yes, what are their new features. Calculations in a deformed single particle 
model [199] have shown that the valence particles in specihc orbitals with low projection 
of the angular momentum on the symmetry axis can give rise to halo structures in the 
limit of weak binding. The deformation of the halo is in this case solely determined by the 
intrinsic structure of the weakly bound orbitals. Indeed, halos in deformed nuclei were 
investigated in several mean held calculations in the past [200, 201, 133]. In Ref. [202], 
it has been concluded that in the neutron orbitals of an axially deformed Woods-Saxon 
potential the lowest-/ component becomes dominant at large distances from the origin 
and therefore all 12^ = 1/2+ levels do not contribute to deformation for binding energies 
close to zero. Such arguments raise doubt about the existence of deformed halos. In 
addition, a three-body model study [203] suggested that it is unlikely to hnd halos in 
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deformed drip line nuclei because the correlations between the nucleons and those due 
to static or dynamic deformations of the core inhibit the formation of halos. 

Therefore a model which provides an adequate description of halos in deformed 
nuclei must include in a self-consistent way the continuum, deformation effects, large 
spatial distributions, and the coupling among all these features. In addition it should 
be free of adjustable parameters in order to make predictions reliable. 

Over the past years, lots of efforts have been made to develop a deformed 
relativistic Hartree (RH) theory [135] and a deformed relativistic Hartree-Bogoliubov 
theory in continuum (DRHBc theory) [108]. As a hrst application, halo phenomena in 
deformed nuclei have been investigated within the DRHBc theory [30]. The detailed 
theoretical framework is available in Ref. [31]. The DRHBc theory is useful to 
answer questions concerning exotic nuclear phenomena in nuclei near the drip line 
[200, 199, 204, 140, 203, 201], such as whether there are deformed halos or not and 
what new features can be expected in deformed exotic nuclei. In Ref. [205], the DRHBc 
theory is extended to incorporate the blocking effect due to an odd nucleon. In such 
a way, pairing correlations, continuum, deformation, blocking effects, and extended 
spatial density distributions in exotic odd-A or odd-odd nuclei can be taken into account 
microscopically and self-consistently. In recent years, the CDFT with the density- 
dependent meson-nucleon couplings have attracted more and more attention owing to 
improved descriptions of the equation of state at high density, the asymmetric nuclear 
matter and the isovector properties of nuclei far from stability. Therefore, it’s necessary 
to develop a density-dependent deformed relativistic Hartree-Bogoliubov (DDDRHB) 
theory in continuum. The challenge is the treatment on the density-dependent couplings, 
meson helds, and potentials in axially deformed system with partial wave method. 
Quite recently, the DDDRHB theory in continuum was developed in Ref. [109]. In 
this Section, the halo in deformed exotic nucleus and its mechanism as well as possible 
shape decoupling will be briefly reviewed. 

6.1. Neutron separation energy and nuclear size 

Figure 18 shows the neutron Fermi energy An and two neutron separation energy S' 2 n of 
Mg isotopes calculated in the DRHBc theory with the parameter sets NL3 and PKl. 
Available data of S' 2 n taken from [183] are also included for comparison. Except the 
different prediction of the drip line nucleus, the neutron Fermi surfaces and two neutron 
separation energies from both parameter sets are very similar. The calculated two 
neutron separation energies S2n of Mg isotopes agree reasonably well with the available 
experimental values except for ^^Mg. The large discrepancy in ^^Mg is connected to the 
shape and the shell structure at A^ = 20. For details, see Ref. [31]. 

Experimentally ^‘^Mg has been observed [206]. Theoretically there are different 
predictions on the last bound nucleus in Mg isotopes, e.g., '^^Mg from the hnite range 
droplet model [207], '^'^Mg from a macroscopic-microscopic model [208], a RMF model 
with the parameter set NLSH [209], and a Skyrme HFB model solved in a 3-dimensional 
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Figure 18 . (Color online) The neutron Fermi energy An (left panel) and two neutron 
separation energy S' 2 n (right panel) of Mg isotopes calculated by DRHBc theory with 
the parameter sets NL3 and PKl. The data of 5'2n (labeled as “Exp.”) are taken 
from [183], Taken from Ref. [31]. 


Cartesian mesh [210], and from a Skyrme HFB model solved in a transformed 

harmonic oscillator basis [64] and the HFB21 mass table [211]. Therefore the prediction 
of the nentron drip line nnclens in Mg isotopes is model and parametrization dependent. 
In the DRHBc calculations with the parameter set NL3, ^®Mg is the last nucleus of which 
the neutron Fermi surface is negative and the two neutron separation energy is positive. 
However, with the parameter set PKl, ^^Mg is predicted to be the last nucleus within 
the neutron drip line. 

In Fig. 19, the rms radii for Mg isotopes are plotted as functions of the neutron 
number, including the neutron radii proton radii i?p, matter radii the 
curve, and the available data for Rt [212, 213]. The calculated matter radius follows 
roughly the curve up to H = 34. From ^®Mg on, the matter radius lies much high 

above the tqA^^^ curve. This indicates some exotic structure in these nuclei. 

6.2. Density distributions and nuclear shape 

Figure 20 shows neutron density prohles of even-even Mg isotopes with H > 28 
calculated in DRHBc theory with the parameter set PKl [31]. The spherical component 
of the neutron density distribution is represented by Pn,A=o(’")- The density distributions 
along the symmetry axis z {6 = 0°) and perpendicular to the symmetry axis z {6 = 90°) 
are respectively represented by Pn{z,r± = 0) and pa{z = 0,rj_) with r± = R y"^. 

With increasing A, the spherical component of the neutron density distribution 
Pn,A=o(^) changes quickly at ^"^Mg and abruptly at ^^Mg. The density distribution along 
the symmetry axis Pni^Nx = 0) changes abruptly from ^^Mg to ^'^Mg. This can be 
understood easily from the shape change in these two nuclei: ^^Mg is spherical but ^'^Mg 
is prolate which means that it is elongated in the .2 axis. In the direction perpendicular 
to the symmetry axis, the neutron density pn(^ = 0,rj_) of ^^Mg extends very far away 
from the center of the nucleus and a long tail emerges, revealing the formation of a halo. 
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Figure 19 . (Color online) The root mean square (rms) radii for Mg isotopes calculated 
in DRHBc theory plotted as functions of the neutron number. The neutron radii 
i?n, proton radii i?p, matter radii i?t, the curve, and the available data for 

i?t [212, 213] are given. The neutron radius i?nj proton radius i?p and total radius Rt 
of Mg isotopes are calculated with PKl. The curve is included to guide the eye. 

The data for matter radii labeled by “Exp.” are taken from Refs. [212, 213]. Taken 
from Ref. [31]. 


By comparing = 0) and p-a_{z = 0,rx) for ^^Mg, it was found that in the 

tail part, the neutron density distributes more along the direction perpendicular to the 
symmetry axis. Since this nucleus as a whole is prolate, it indicates that the tail part 
has a different shape from the nucleus itself. This is the shape decoupling which was 
hrst predicted in Ref. [30]. 

The comparison of quadrupole deformation f3 between the deformed RHB 
calculation and the experiment is given in Fig. 21. The experimental value of (3 is 
extracted from the measured B{E2 : O]*' 2f) and only the absolute value could be 

obtained [214]. The results of NL3 are almost the same as the results of PKl except for 
^°Mg. R is slightly oblate with NL3 but it has a prolate shape with PKl. For ^^Mg, the 
gap between the neutron levels 1 ^ 3/2 and I/ 7/2 is almost 7 MeV which results in a strong 
closed shell at = 20. Therefore the deformed RHB calculations with both parameter 
sets predict spherical shapes for this nucleus. This results in the large discrepancy in 
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Figure 20 . (Color online) Neutron density profiles of even-even Mg isotopes with 
A > 28 calculated in DRHBc theory with the parameter set PKl. The upper panel 
represents the spherical component of the neutron density distribution Pn,x=o{f)- The 
middle panel refers to the density distribution Pn(z, c_l = 0) along the symmetry axis 
z { 0 = 0°) with r± = +y^ and the lower panel refers to pn{z = 0,r_L) that 

perpendicular to the symmetry axis z {0 = 90°). Taken from Ref. [31]. 
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Figure 21 . (Color online) The quadruple deformation parameter /3 of Mg isotopes 
calculated with the parameter sets NL3 and PKl. The data labeled with “Exp.” are 
taken from [214]. Taken from Ref. [31]. 


S' 2 n for as is seen in Fig. 18. Other mean field models also predict spherical or 

almost spherical shapes for ^^Mg [215, 209, 210, 216, 217]. However, this nucleus is well 
deformed as is indicated by the large B{E2 : 0]'' —?■ 2]*') value [218]. The same is true 
for some neighboring nuclei with iV ~ 20. These nuclei form the so called “island of 
inversion” [219, 220, 221, 222] which is related to the quenching of the iV = 20 shell 
closure. In order to describe the shell quenching effect at iV = 20, one has to go beyond 
the mean held [223, 217]. 

6.3. Deformed halo and shape decoupling 

The deformed relativistic Hartree-Bogoliubov theory in continuum with the parameter 
set PKl predicted the last bound neutron rich nucleus for Mg isotopes as ^^Mg [31]. 
There are both prolate and oblate stable solutions for "^^Mg. The ground state of ^^Mg 
is prolate. With the triaxially deformed relativistic mean held calculations [224], it is 
shown that the obtained oblate minimum is not an isomeric state but a saddle point in 
the potential energy surface. 

It is found that the ground state of ^^Mg is well deformed with fd ~ 0.41 and two 
neutron separation energy S 2 n ~ 0.22 MeV [31]. The density distribution of such weakly 
bound nucleus has a very long tail in the direction perpendicular to the symmetry axis, 
see Fig. 20, which indicates the prolate “^^Mg has an oblate halo. 

The density distribution is divided into contributions of the oblate “halo” and of 
the prolate “core” in Fig. 22. The density distribution of this weakly bound nucleus has 
a very long tail in the direction perpendicular to the symmetry axis which indicates the 
prolate nucleus “^^Mg has an oblate halo and there is a decoupling between the shapes 
of the core and the halo. It should be noted that in the calculations with NL3, “^^Mg is 
the last nucleus [30]. 

The single particle spectrum around the Fermi level for the ground state of ^^Mg 
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Figure 22 . (Color online) Density distributions of the ground state of with the 
z axis as the symmetry axis: (a) the neutron halo, and (b) the neutron core. Taken 
from Ref. [31]. 


is shown in Fig. 23 [31]. The good qnantnm nnmbers of each single particle state are 
also shown. The occnpation probabilities in the canonical basis have BCS-form [188] 
and are given by the length of the horizontal lines in Fig. 23. The levels close to the 
threshold are labeled by the nnmber i according to their energies, and their conserved 
qnantnm nnmber VC as well as the main spherical components are given at the right 
hand side. The nentron Fermi level is within the pf shell and most of the single particle 
levels have negative parities. Since the chemical potential A„ is close to the continnnm, 
orbitals above the threshold have noticeable occnpations dne to the pairing correlations. 
The single nentron levels of “^^Mg can be divided into two parts, the deeply bonnd levels 
(^can < —2 MeV) corresponding to the “core”, and the remaining weakly bound levels 
close to the threshold (^can > —0.3 MeV) and in the continuum corresponding to the 
“halo”. 

As discussed in Refs. [30, 31, 225], the shape of the halo originates from the intrinsic 
structure of the weakly bound or continuum orbitals. By examining the neutron density 
distribution, it was found that for the ground state of '^^Mg, the halo is mainly formed 
by level 4 and level 5. A detailed analysis of the intrinsic properties revealed that in 
'^^Mg the shape of the halo is oblate and decouples from the prolate core. 

We note that this kind of shape decoupling was later also found in non-relativistic 
HFB calculations [226]. An exotic “egg”-hke structure consisting of a spherical core plus 
a prolate halo was predicted in ^®Ne, in which the near-threshold nonresonant continuum 
plays an essential role. 

6.4- Recent progress 

For odd particle system, the halo could depend strongly on the interplay among the 
odd-even effects, continuum and pairing effects, deformation effects, etc. Combining 
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Figure 23 . (Color online) Single neutron levels of ground state of ^^Mg in the 
canonical basis as a function of the occupation probability The blue dashed line 
corresponds to the BCS-formula with an average pairing gap. Taken from Ref. [31]. 


the necessity for the improved descriptions of the equation of state at high density, 
the asymmetric nuclear matter and the isovector properties of nuclei far from stability 
with the density-dependent meson-nucleon couplings [227, 228, 229, 230, 138, 182], it 
is necessary to extend the DRHBc theory with the density-dependent meson-nucleon 
couplings. Here we will briefly review these progress, i.e., the exotic nuclear structure in 
unstable odd-H or odd-odd nuclei [205] and the density-dependent deformed relativistic 
Hartree-Bogoliubov (DDDRHB) theory in continuum [109]. 

In order to describe the exotic nuclear structure in unstable odd-H or odd-odd 
nuclei, the blocking effect of one or several nucleons has to be taken into account. In 
Ref. [205], the DRHBc theory [108, 30, 31] is extended to incorporate the blocking effect 
due to odd particle system. In such a way, pairing correlations, continuum, deformation, 
blocking effects, and extended spatial density distributions in exotic odd-H or odd-odd 
nuclei can be taken into account microscopically and self-consistently. 

For an odd nucleon system with the Fth level blocked, the time reversal symmetry 
is violated and there appear currents in the system. In principle, one has to diagonalize 
the RHB equation twice. Using the equal hlling approximation which is usually 
made [231, 232], the currents are neglected and the two conhgurations of a particle in the 
time-reversal partner space are averaged in a statistical manner [205] . The corresponding 
currents cancel each other and in this way, in each step of the iteration meson helds 
with the time reversal symmetry is obtained. In practice, the density matrix p and the 
abnormal density k in two subspaces are averaged and replaced by 

p' = i {UkMl - KK) . 


( 69 ) 
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Figure 24 . (Color online) Neutron single particle levels in the canonical basis and 
their occupation probability for from the DRHBc and RCHB calculations. 
Taken from Ref. [205]. 


- 5 + v-jjl ), ( 70 ) 

where Vk^ and f/^ are column vectors in the matrices V and U corresponding to the 
blocked level. 

Using the density functional PKl [138] and a zero range density dependent pairing 
force [39, 31], the deformed RHB equations are solved in a spherical Dirac Woods- 
Saxon basis [65] and the deformed potentials and densities are expanded in terms of the 
Legendre polynomials Px{cos9) (for details see Refs. [108, 30, 31]). In order to check 
the accuracy of the DRHBc code with the blocking effect, the spherical nucleus is 
investigated and the results of the DRHBc code allowing only the spherical components 
of the helds, i.e., A = 0, are compared with those obtained in RCHB theory [39]. It is 
found that both calculations agree well with each other with high precision [205]. 

In Fig. 24, the neutron single particle levels in the canonical basis for is shown 
in comparison with the RCHB results in order to understand the details and validity 
of the DRHBc calculation. The length of each level is proportional to the occupation 
probability. In the RCHB calculation, the blocked orbital is 1 ^ 5 / 2 - In the DRHBc code, 
a spherical solution is enforced by allowing only A = 0 component in the Legendre 
expansion mentioned earlier and therefore the three sublevels of the 1 ^ 5/2 orbital with 
= 1/2+, 3/2+, and 5/2+ are degenerate. The present results in Fig. 24 are obtained 
with the level = 1/2+ blocked in the DRHBc code, but blocking any of these three 
levels gives the same results. From Fig. 24 one hnds nice agreement between the results 
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Figure 25 . (Color online) (a) Neutron density profiles of and in DRHBc 
calculations with the parameter set PKl. (b) Neutron single particle levels in the 
canonical basis and their occupation probability for in DRHBc calculation. 
Taken from Ref. [205]. 


from the DRHBc and RCHB calculations. 

An application of the DRHBc theory for an exotic and deformed nucleus with an 
odd number of particles is given in Ref. [205] for the deformed halo candidate nucleus 
[233] with the parameter set PKl [138]. In contrary to the spherical even-even 
nucleus the DRHBc theory predicts a deformation (3 = 0.25 for the odd-A nucleus 
The neutron root mean square radii of and are 2.56 fm and 2.79 fm, 
respectively. 

The DRHBc calculations for lead to a deformed density distribution with an 
axial symmetry. In Fig. 25, the neutron densities along the symmetry axis Pni'r'z, f'± = 0) 
and perpendicular to the symmetry axis Pni^z = 0, r±) are plotted as dashed and dashed- 
dotted lines. The spherical density for is also included as a reference by a solid line. 
It is interesting to see that, in the direction perpendicular to the symmetry axis, the 
neutron densities are almost the same for and at least when < 6 fm. Along 
the symmetry axis the neutron density of extends much further than that of 
This is partly due to that is prolate and the weakly bound 1/2+ level is occupied. 

The single particle levels of in the canonical basis are plotted in Fig. 25. Since 
it is a deformed nucleus, the lp 3/2 orbital is split into two levels with = 1/2“ and 
= 3/2“, respectively. There is one neutron occupying in = 1/2+ level near the 
threshold. The occupation probability = 0.5, indicating that it is averaged over the 
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two configurations with = ±1/2+. Because of the deformation, this level is a mixture 
of the spherical orbitals 1 ^ 5/2 (62%) and 2si/2 (36%). The weakly-bound feature and 
the relatively large s-wave component of this level results in that the neutron density of 

extends further along the symmetry axis. 

After the comparison between the DRHBc and the spherical RCHB results for 
as well as the application for the deformed neutron-rich nucleus the DRHBc with 
blocking effect proved to be a powerful tool for the study of the odd-A or odd-odd exotic 
nuclei. 

In the development of the DRHBc with the density-dependent meson-nucleon 
couplings [109], the challenge is the treatment of the density-dependent couplings, meson 
helds, and potentials in axially deformed system with partial wave method. 

In Fig. 26, by taking the a meson as an example, the partial waves of 
density-dependent coupling strengths for ^®Mg are shown with g^iPv) = 

XIa 5 'o',a(^)-Pa(cos 6 '). It can be seen that the major component is that with A = 0, 
which is more than one order of magnitude larger than the others. The amplitudes of 
gcr,x decrease quickly with increasing A and become negligible when A > 8 . This means 
that a cutoff Xg^ranx = 10 is enough to provide accurate results. 



0 2 4 6 8 

r(fm) 

Figure 26 . (Color online) Decompositions of coupling strength g^-.A with A = 0 (a) 
and A = 2,4, 6,8,10 (b) for ^®Mg calculated with density-dependent DRHBc theory. 
Taken from Ref. [109]. 
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By comparing the calculated bulk and single particle properties of neutron-rich 
nucleus ^®Mg in spherical case with Ap^max = 0 with those obtained by RCHB [39], 
the density-dependent DRHBc theory is verihed [109]. It therefore opens the door to 
investigate the deformed exotic nuclei with effective density-dependent interactions. 

7. Summary and perspectives 

In this Topical Review, we have reviewed the progresses on the covariant density 
functional theory in continuum for the neutron halo phenomenon in both spherical 
and deformed nuclei. 

In the mean held study of nuclear halos, a proper description of nuclear pairing and 
the continuum contribution is very important. To include properly the contribution of 
continuum states, the Bogoliubov transformation has been justihed to be very useful. 
The relativistic continuum Hartree-Bogoliubov (RCHB) theory has been developed in 
which the coupled differential equations are solved in coordinate (r) space and the mean 
held ehects of the coupling to the continuum can be fully taken into account. 

From the RCHB theory, a self-consistent description of the neutron halo in ^^Li 
has been achieved. The calculated density distribution agrees well with the experiment. 
It has further been shown that the halo in ^^Li is formed by Cooper-pairs scattered 
in the two levels lpi /2 and 2si/2 with the former below the Fermi surface and the 
latter in the continuum. In the RCHB theory, a halo appears if there are several single 
particle valence orbitals with small orbital angular momenta and correspondingly small 
centrifugal barriers close to the continuum limit. 

For medium and heavy nuclei, the RCHB theory has treated the influence of 
pairing correlations and many-body effects properly and predicted the existence of more 
neutrons in the halo. Giant neutron halos, consisting of up to six halo neutrons, have 
been predicted in Zr and Ca isotopes close to the neutron drip line. When approaching 
the neutron drip line, the two-neutron separation energies for several nuclei in Ca and 
Zr isotopes are extremely close to zero, which means that the valence neutrons gradually 
occupy the loosely bound states and the continuum. 

The phenomena of halo and giant halo of neutron-rich even-Ca isotopes were 
investigated and compared in the framework of the RCHB and non-relativistic Skyrme 
HFB calculations. It has been found that although halo phenomena exist for Ca isotopes 
near neutron drip line in both calculations, the halo in the Skyrme HFB calculations 
starts at a more neutron-rich nucleus than that in the RCHB calculations, and the RCHB 
calculations have larger neutron root-mean-square radii systematically in iV > 40 than 
the Skyrme HFB calculations. The former difference comes from the difference in shell 
structure. The reasons for the latter can be partly explained by the neutron 3si/2 orbital, 
which causes more than 50% of the difference for neutron radii at ®®Ca. 

The relativistic Hartree-Fock-Bogoliubov theory in continuum (RHFBc theory) 
with density-dependent meson-nucleon couplings has been developed. With this theory, 
the influence of exchange terms on the formation of neutron halos has been studied. 
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Giant halos as well as ordinary ones were found in Ce isotopes close to the neutron drip 
line. It was demostrated that the shell structure evolution, when approaching the drip 
line, is important for the appearance of neutron halos. 

The RCHB theory has also been extended to including hyperons and even Ca 
isotopes ranging from the proton drip line to the neutron drip line have been studied 
systematically. Compared with ordinary nuclei, the additional one or two A-hyperons 
lower the Fermi level and the last bound hypernuclei have two more neutrons than the 
corresponding ordinary nuclei. Giant halo phenomena due to the pairing correlation and 
the contribution from the continuum have been predicted in Ca hypernuclei, similar to 
those appearing in ordinary Ca isotopes. 

A deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc theory) 
has been developed aiming at a proper description of exotic nuclei, particularly those 
with a large spatial extension. In order to give an adequate consideration of both the 
contribution of the continuum and the large spatial distribution in exotic nuclei, the 
deformed RHB equations were solved in a Woods-Saxon (WS) basis in which the radial 
wave functions have a proper asymptotic behavior at large distance from the nuclear 
center. 

Halo phenomena in deformed nuclei were investigated within the DRHBc theory. 
These weakly bound quantum systems present interesting examples for the study of 
the interdependence between the deformation of the core and the particles in the halo. 
The Mg and Ne isotopes have been studied in details. It was found that in 
there appears a decoupling of the halo orbitals from the deformation of the core: The 
core is prolate, but the halo has a slightly oblate shape. The generic conditions for the 
occurrence of this decoupling effects have been given: the existence and the deformation 
of a possible neutron halo depends essentially on the quantum numbers of the main 
components of the single particle orbitals in the vicinity of the Fermi surface. 

Dehnitely there exist lots of open questions on the study of nuclear halos with the 
covariant density functional theory. We list several of them below. 

In the RCHB, RHFBc and DRHBc theories, the box boundary condition has been 
applied and the continuum is discretized. Therefore one cannot get the width for 
resonances. With the Green’s function method, the correct asymptotic behavior can be 
given on the wave function for the continuum states. The Green’s function method has 
been implemented in the non-relativistic HFB theory [55]. The effort to implement the 
Green’s function approach in the RHB theory is undertaken and promising preliminary 
results have been obtained [58]. 

In the systematic study of giant halos in severl isotope chains, it has been found 
that the appearance of giant halos and the contribution of continuum are very important 
in determining the boundary of nuclear chart. Furthermore, the deformation effects will 
also influence very much the position of drip lines. More systematic study is now carried 
out to reveal the global impact of the continuum and deformation effects on the nuclear 
boundary. 

If the core and the halo of a nucleus have different shapes, many interesting 
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questions arise: For example, What will happen when the nucleus rotates? What 
role could the shape fluctuation play? To answer these questions, future works are 
welcome to implement the angular momentum projection in the DRHBc theory and to 
use the generator coordinate method (GCM) or to solve Bohr Hamiltonian based on the 
potential energy surface obtained from constraint DRHBc calculations. 

The possible experimental signals of shape decoupling in deformed halo nuclei 
are certainly interesting topics. The halo feature is connected with relatively large 
cross sections and narrow longitudinal momentum distributions in knockout reactions. 
The decoupling between the deformations of the core and the halo may manifest itself 
by some new experimental observables, e.g., the double-hump shape of longitudinal 
momentum distribution in single particle removal reactions and new dipole modes, etc. 
In particular, a combination of the experimental method proposed in Ref. [234] and the 
theoretical approach developed in Ref. [235] would be useful in the study of longitudinal 
momentum distribution in single particle removal reactions with deformed halo nuclei as 
projectiles. The shape decoupling effects may also has some influence on the sub-barrier 
capture process in heavy ion collisions [236]. 
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